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ABSTRACT 


OPTIMAL  SHIP  POSITIONS  FOR 
NATAL  BATTLE  GROUP  DEFENSE  PROBLEMS 


Robin  Cynthia  Magonet-Neray 


Submitted  to  the  Department  of  Electrical  Engineering 
on  March  16,  1983  in  partial  fulfillment  of  the 
,  requirements  for  the  Degree  of  Master  of  Science  in 

\  .  Operations  Research 

T?his  report  presents  an  optimization  model  to  maximize  the  survival  probabili¬ 
ty  of  a  carrier  operating  in  a  Battle  Group  (BG)  environment  given  AAV 
(anti-air  warfare)  and  ASM  (anti-sub  warfare)  resources.  The  model  is  a  stat¬ 
ic,  probabilistic,  two-dimensional  representation. 

/'""W :*ev*iep<j^everal  probabilistic  and  parametric  relationships,  including: 

(ly  the  probability  that  a  ship  hills  a  target 
_ [2^  the  optimum  intercept  distance 
-v3>'the  effect  of  the  ships'  distance  to  the  carrier 
u.  Jtf  the  probability  distributions  of  the  targets, 

"'fcvjA  survival  probability  of  the  carrier, 

Ci  The  solution  to  the  problem  is  the  optimum  location  of  the  AAW  and  ASW  ships 
with  respect  to  the  carrier;  those  locations  that  maximize  the  probability  of 
survival  of  the  carrier,  from  the  AAW  and  ASW  threats.?-"' 

We  perform  numerical  experiments  when  both  one  and  many  ships  are  defending 
the  carrier.  We  study  the  quantitative  impact  on  the  survival  probability  of 
the  carrier  and  on  the  optimum  location  of  the  ship(s)  when  we  do  parametric 
analyses  that  reflect  changes  in 

(a)  threat  sector  locations 

(b)  target  density  in  each  sector 

(c)  AAW  and  ASW  defense  capabilities  on  a  per  ship  basis 

(d)  number  of  defending  ships 
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CHAPTER  1  INTRODUCTION 


1.1  MOTIVATION  AND  PURPOSE 


With  the  rising  cost  and  sophistication  of  aircraft  carriers,  it  is  increas¬ 
ingly  more  important  to  find  better  ways  to  protect  them. 


Typically,  a  carrier  is  defended  against  air  and  submarine  threats  by  ships 
with  simultaneous  AAW  (anti-air  warfare) ,  ASW  (anti-submarine  warfare)  and 
ASUW  (anti-surface  warfare)  capabilities.  These  ships  are  located  in  some 
"neighbourhood"  of  the  carrier.  With  the  aid  of  advanced  radar  and  sonar, 
these  ships  are  required  to  detect  and  kill  enemy  air,  surface  and  sub 
targets.  Clearly,  the  locations  of  these  ships  relative  to  the  carrier  deter¬ 
mines  their  effectiveness  against  enemy  threats. 


The  purpose  of  this  report  is  to  develop  a  preliminary  mathematical  model  to 
study  such  problems.  For  the  sake  of  simplicity,  we  focus  attention  only  on 
the  AAW  and  ASW  threats.  The  model  is  a  static,  probabilistic  and 
two-dimensional  representation.  Its  goal  is  to  maximize  the  probability  of 
survival  of  the  carrier  given  air  and  sub  threats,  by  finding  the  optimum 
locations  of  the  defending  ships,  relative  to  the  carrier.  We  test  the  per¬ 
formance  of  the  model  through  comprehensive  numerical  studies. 


Chapter  l  Introduction 


We  wish  to  stress  at  the  outset  that  the  mathematical  model  as  presented  here 
cannot  be  used  for  realistic  ship  location  decision-aids  for  an  actual  battle 
group.  The  models  that  we  have  used  are  too  crude  to  properly  capture  the  AAW 
and  ASW  surveillance  and  weapons  effectiveness  of  actual  platforms 
(destroyers,  frigates,  cruisers,  etc.).  Also,  the  numerical  values  that  have 
been  used  in  the  models,  do  not  represent  any  real  platform  parameters  (which 
are  classified) .  The  real  intent  of  this  research  was  to  under stand  the  com¬ 
plexities  of  the  optimization  problem  involved  if  a  more  realistic  model  were 
to  be  used. 

The  Battle  Group  location  problem  is  a  very  important  problem  which  is  amena¬ 
ble  to  quantitative  analysis  and  optimization,  Athans  (l) ,  (2).  However,  this 
problem  has  not  been  analyzed  in  the  unclassified  literature,  except  for  the 
very  recent  study  reported  by  Castanon  et  al  (3),  which  is  in  the  same  spirit 
as  the  research  reported  herein.  However,  the  quantitative  defense  effective¬ 
ness  models  in  (3)  are  quite  different. 

In  spite  of  the  limitations  of  these  models,  they  offer  valuable  guidance  on 
how  to  exploit  the  multiwarfare  capabilities  of  modern  naval  platforms  and  to 
use  them  in  a  coordinated  fashion  in  the  defense  of  a  modern  Battle  Group.  We 
quote  below  a  paragraph  of  the  paper  (4)  by  R.F.  Schoultz  (RADM,  USN)  that 


"  The  situation  for  surface  combatants  is  perhaps  the  most 

illustrative.  With  the  shrinking  of  the  force  levels  there  is  a 
need  to  share  platforms  among  the  missions  of  AAW,  asw,  and  ASUW. 
Depending  upon  the  scenario,  these  forces  must  be  allocated  to 
maximize  their  effectiveness  against  the  predominate  threat 
concentration  at  the  given  time...." 

1.2  SCOPE  AND  ORGANIZATION 

In  Chapter  2  we  state  and  formulate  the  problem  which  is  the  focus  of  this 
report.  The  problem  is  to  maximize  the  probability  of  survival  of  a  carrier 
given  air  and  sub  threats,  under  different  scenarios.  The  solution  is  the 
optimum  location  of  the  defending  ships  relative  to  the  carrier;  those 
locations  that  maximize  the  survival  probability  of  of  the  carrier.  We  build 
the  model;  we  present  the  principles,  define  the  parameters  and  develop  the 
algebraic  and  probabilistic  relationships.  We  address  the  analytic  limita¬ 
tions  of  the  model. 

The  purpose  of  Chapter  3  is  to  study  the  behaviour  of  the  model  when  a  single 
ship  is  defending  the  carrier.  Our  interest  lies  in  discovering  the  effect  on 
the  survival  probability  of  the  carrier  and  on  the  optimum  location  of  the 
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ship  when  we  do  a  parametric  analysis.  We  examine  the  consequences  of  changing 
the  sector  and  increasing  the  number  of  targets.  We  explore  the  existence  of 
local  maxima. 

In  Chapter  4  we  further  test  the  model  by  considering  scenarios  in  which  two 
and  three  ships  are  protecting  the  carrier.  Our  experimental  studies  focus  on 
determining  how  multiple  ships  work  in  conjunction  with  one  another  to  best 
defend  the  carrier.  We  investigate  the  conditions  under  which  symmetry  exists 
and  consider  the  result  of  pure  (AAW  or  ASW  threats,  but  not  both)  threats.  We 
study  the  effect  of  multiple  targets,  probe  the  limiting  cases,  and  examine 
the  consequences  of  changing  the  threat  sectors. 

Chapter  5  proposes  natural  extensions  to  the  model  and  is  aimed  at  stimulating 
the  mind  of  the  reader  to  invent  his  own  projects. 

We  summarize  conclusions  in  Chapter  6. 
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CHAPTER  2  PROBLEM  STATEMENT  AND  MATHEMATICAL  FORMULATION 


2.0  SUMMARY 


The  purpose  of  this  chapter  is  to  define  the  optimization  problem  which  is  the 
focus  of  this  report,  and  to  develop  a  mathematical  model  which  captures  the 
essence  of  the  problem.  The  problem  is  to  maximize  the  probability  of  survival 
of  a  carrier  given  air  and  sub  threats,  under  various  scenarios.  The  solution 
is  the  optimum  locations  of  the  ships  defending  the  carrier. 

The  model  is  a  static,  probabilistic  and  two-dimensional  representation.  The 
probability  that  a  ship  intercepts  a  target,  the  optimum  intercept  distance, 
the  effect  of  a  ship's  distance  to  the  carrier  and  the  probability  that  the 
carrier  survives,  are  all  developed.  The  limitations  of  the  model  are 
addressed. 


2.1  PROBLEM  STATEMENT 


The  purpose  of  the  problem  is  to  maximize  the  probability  of  survival  of  a 
carrier  given  air  and  sub  threats,  given  simple  models  of  the  survival  proba¬ 
bility. 
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There  are  AAW  (anti-air  warfare)  and  ASW  (anti- sub  warfare)  ships  available  to 
defend  the  carrier.  Missiles  (air  threats)  and  submarines  (sub  threats)  are 
approaching  the  carrier  uniformly  over  well-defined  sectors,  which  are  known 
with  certainty.  Each  of  the  AAW  and  the  ASW  ships  (cruisers  and  destroyers 
respectively)  have  both  anti-air  and  anti-sub  capabilities.  That  is,  an  AAW 
ship  has  superior  anti-air  warfare,  but, is  still  capable  and  responsible  for 
intercepting  sub  targets.  Similarly  for  ASW  ships. 

The  solution  to  the  problem  is  the  optimum  location  of  the  ships  with  respect 
to  the  carrier  -  that  location  which  maximizes  the  probability  of  survival  of 
the  carrier. 

We  require  that  ships  not  be  positioned  at  the  origin  (the  assumed  location  of 
the  carrier)  since  it  is  assumed  that  a  certain  distance  is  necessary  to  suc¬ 
cessfully  detect  the  targets  (surveillance).  This  constraint  is  incorporated 
into  the  mathematical  model.  Figure  2.0  illustrates  this  scenario. 
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2.2  THE  MODEL  -  A  MATHEMATICAL  FORMULATION 


The  problem  is  formulated  in  2-dimensions.  The  carrier  is  positioned  at  the 
origin.  Each  ship  i  has  Cartesian  coordinates  (x(*  ,y(-  )  which  represent  its 
position,  relative  to  the  carrier.  We  assume  that  threats  are  approaching  the 
carrier  at  an  angle  9  which  is  uniformly  distributed  over  a  known  sector. 


2.2.1  OPTIMUM  INTERCEPTION  DISTANCE 


The  ship  will  attempt  to  intercept  a  single  target  when  the  distance  between 
them  is  at  a  minimum.  This  is  the  optimum  interception  distance,  and  we  now 
establish  how  it  is  determined. 

Let  the  (fixed)  position  of  ship  i  be  (xj ,y* ) .  The  target  is  approaching  the 
carrier  at  an  angle  9,  and  at  any  given  instant  has  as  its  location  coordi¬ 
nates  (Rcosfl ,Rsin0) ,  where  R  is  the  current  radius*  (distance)  of  the  target  to 
the  carrier.  Let  d;  be  the  distance  of  ship  i  to  the  target.  We  wish  to  find 
d*  ,  the  optimum  interception  distance. 

x 

The  Euclidean  distance  dj  (see  Figure  2.1)  is: 

X  2.  Xj 

dj  =  (x;  -  Rcosfl)  +  (yj  -  Rsinfl)  (2.1) 

The  optimum  distance  d*  occurs  at  =  0,  which  yields 
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Substituting  the  value  of  R  in  Eq'n  (2.2)  into  Eq'n  (2.1),  we  obtain 

d*  =  (XjSintf)  -  (y,-cos0)  (2.3) 

which  is  the  optimum  distance  for  ship  i  to  intercept  a  target. 

DEFINITION  2.1:  d* 

d*j  is  the  INTERCEPT  DISTANCE  and  is  given  by 
d*,-  =  (xj  sin0)  -  (yj  cos0) 

Figure  2.2  shows  d*  . 


2.2.2  SOME  PROBABILISTIC  RELATIONSHIPS 


Each  cruiser  and  destroyer  has  characteristic  abilities  to  kill  air  and  sub 
targets.  These  air  and  sub  capabilities  are  described  in  terms  of  probabili¬ 
ties  of  kill,  and  we  assume  that  they  are  a  function  of  both  the 
ship-to-target  intercept  distance,  d*.  ,  and  the  distance  r.  of  the  ship  to 
the  carrier,  r;  =4xj*"+  yj*. 
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2. 2. 2.1  THE  EFFECT  OF  A  SHIP'S  DISTANCE  TO  THE  CARRIER 

Since  the  carrier  is  fixed  at  the  origin,  all  targets  of  interest  pass  through 
the  origin. 

In  Section  2.2.1,  we  saw  that  the  ship  will  attempt  to  kill  the  target  when 
the  distance  between  them,  d,-  ,  is  at  a  minimum.  If  minimizing  this  distance 
were  the  only  objective,  then  the  optimum  solution  would  suggest  placing  all 
the  ships  at  the  origin,  since  then  d?  of  Eq'n  (2.3)  would ^equal  zero.  This  is 
clearly  infeasible  since  the  ships  cannot  be  on  top  of  the  carrier.  It  is  also 
undesirable  since  the  ships  must  maintain  a  certain  distance  from  the  carrier 
in  order  to  detect  the  targets  at  a  sufficiently  large  distance  from  the  car¬ 
rier. 

To  incorporate  this  constraint  into  the  problem,  we  force  the  probability  of 
intercepting  a  target  to  be  a  function  of  the  ship's  distance  to  the  origin. 
We  call  this  function  Pr  .  The  probability  of  ship  i  killing  a  target 
increases  as  its  distance  to  the  carrier  increases.  This  prevents  the  ship 
from  having  its  final  position  at  the  origin. 

This  assumption,  that  the  further  the  ship  is  from  the  carrier  the  better  is 
the  chance  that  it  will  kill  a  target  (for  the  same  intercept  distance) 
appears  to  make  sense.  ASW  ships,  employing  passive  sonar,  cannot  be  located 
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near  the  carrier  because  of  noise  interference.  For  AAW  ships  we  may  have  sim¬ 
ilar  electromagnetic  interference  phenomena;  also,  we  want  to  avoid  the  battle 
front  being  too  close  to  the  carrier. 

OEFIMTIOH  2.2:  /> 

'a; 

P,  .  ,  is  defined  as: 

rAi 

-Ofc  r,-  ) 

\ •  1  -  e  ,  k4l>  0  (2.4) 


where , 

r-  =  y  *-  is  the  (Euclidean)  distance  of  ship  i  to  the  carrier 

at  the  origin,  and 

Tk.a.  is  a  parameter  associated  with  each  cruiser  and  destroyer. 
Figure  2.3  shows  the  shape  of  Pr^  . 


The  value  of  selected  affects  the  optimal  position  of  the  ships  -  a 
smaller  k^j  value  results  in  a  ship  being  farther  from  the  carrier,  k^.  is 
chosen  so  as  to  place  the  ships'  distances  from  the  carrier  within  a  reason¬ 
able  range  (i.e.  20  -  100  miles).  In  addition,  the  k^  value  for  destroyers 
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Pr^.  tends  to  1  as  r-  ,  the  distance 

'A  * 

infinity.  It's  smallest  value  is  "0" 
ship  i  to  the  carrier  is  zero.  This  is 
the  ships  away  from  the  origin. 

2. 2. 2. 2  PROBABILITY  OF  INTERCEPTION 

Each  ship  i  has  two  functions  which  fully  describe  its  ability  to  to  intercept 
targets.  PAj  of  Eq'n  (2.5)  reflects  ship  i’s  ability  to  intercept  AAW  tar¬ 
gets.  An  analogous  function,  P^ j  ,  reflects  ship  i’s  ASW  intercept 
capability. 

P-.  and  P-.  are  probabilites  since  their  values  lie  in  (0,1).  However,  they 
are  not  probability  distributions  since  they  do  not  sum  to  unity  as  they  are 
not  scaled  accordingly  . 

Every  ship  i  has  its  own  range  for  killing  both  AAW  and  ASW  targets.  This 
suggests  the  following  definitions. 


from  the  ship  to  the  carrier,  tends  to 
and  is  attained  when  the  distance  from 


how  we  use  Pr^.  in  our  model  to  force 
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0£ FI  MIT  I  OH  2.3:  IA;  ,  I5{ 


£a;  is  called  the  AAW  RANGE  FACTOR  of  ship  i  and  represents  ship  i's 
ability  to  hill  AAW  targets. 

I- •  is  called  the  ASV  RANGE  FACTOR  of  ship  i  and  represents  ship  i's 
ability  to  hill  Asw  targets. 


DEFINITION  2.4 :  PAl 


PA>t  ,  the  PROBABILITY  (SHIP  I  KILLS  AN  AAW  TARGET),  is 


(-d*.  /X^  ) 


where , 


P,  *  i  -  e 

rA\ 


-(K*r  <) 


,  y  o 


d*  *  (x,-sinfl  -  yjcosfl) 


y is  the  range  AAW  range  factor. 


(2.5) 


(2.4) 


(2.3) 


Figure  2.4  shows  PA>.  . 

We  see  that  P^  ,  the  probability  that  ship  i  kills  an  AAW  target,  decreases 
exponentially  with  the  distance  d»-  .  In  addition,  the  probability  of  kill 
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increases  as  r-  ,  the  distance  from  the  ship  to  the  carrier  increases.  A  larg¬ 
er  value  is  favourable  and  yields  a  larger  P^.  . 

Similarly,  for  every  ship  i  we  have  a  function  which  completely  describes  its 
ability  to  intercept  ASW  targets.  We  call  this  function  P-.  . 

^>i 

DEFINITION  2.5:  Psi 


P..  ,  the  PROBABILITY (SHIP  I  KILLS  AH  ASW  TARGET),  is 


(-<!**■ /I5 .  ) 


where , 


) 


d«.  =  (xt-sin0  -  y*  cos#) 


and, 

I..  is  the  ASW  range  factor. 


P^.  has  exactly  the  same  properties  as 

Si 


(2.6) 


(2.4) 


(2.3) 
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2. 2. 2. 3  PROBABILITY  DISTRIBUTION  OF  THE  TARGET 


The  targets  are  approaching  the  carrier  at  an  angle  d  which  we  assume  is  uni¬ 
formly  distributed  over  a  given  sector.  The  sectors  are  described  in  angles. 
Let  f(d)  be  the  probability  density  function  (p.d.f.)  of  a  single  target. 
Then  f(0)  is  given  by: 

f (8)  *  l/(b  -  a)  a  3  8  4  b 

=  0  otherwise 

• 

Figure  2.5  shows  f(0). 


2.2.3  SURVIVAL  PROBABILITY  OF  THE  CARRIER 


We  now  have  all  the  information  required  to  construct  the  probability  function 
associated  with  the  survival  of  the  carrier.  We  initially  consider  the  case  of 
a  single  ship  defending  the  carrier,  and  a  ‘ single  target  attacking  the 
carrier,  and  later  expand  to  include  the  possibility  of  multi-ship  protection 
and  multi-target,  multi-sector  attacks.  We  assume  throughout  that  the  carrier 
has  no  defensive  capabilities  of  its  own. 
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2. 2. 3.1  PERMISSIBLE  SECTORS  AND  LIMITATIONS  OF  THE  MODEL 


In  this  section,  we  present  two  hypothetical  scenarios  based  on  two  different 
sectors.  We  show  how  the  expression  for  the  probability  of  survival  of  the 
carrier  changes  when  different  types  of  sectors  are  used. 

Case  la:  1  Ship ,  1  Target,  1  Sector 

We  formulate  the  problem  #in  terms  of  an  AAW  threat,  i.e.  an  AAW  target  is 
approaching  the  carrier  at  an  angle  6  that  is  uniformly  distributed  over  a 
sector  defined  by  (a,b).  All  calculations  for  an  ASW  target  are  analogous. 
Consider  the  scenario  pictured  in  Figure  2.6. 

For  this  scenario  we  have  that: 

Prob(ship  i  kills  an  AAW  target)  =  P^-  (2.5) 

so  that  , 

Prob(ship  i  kills  an  AAW  target  over  sector  (a,b)) 

=  \b  f(0)  x  PA.  d0  (2.7) 

Jq  a* 

This  implies  the  following  about  the  survival  probability  of  the  carrier. 
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T W<f€T 

tr/v:k^ 


© 


(2.7) 


Prob(carrier  survives)  =  Prob(the  ship  kills  the  AAW  target) 

f(tf)  x  Pa-  dtf 

Figure  2.7  shows  the  area  of  interest  of  the  survival  probability  of  the  car¬ 
rier. 

We  now  establish  the  survival  probability  of  the  carrier  under  a  slightly  dif¬ 
ferent  scenario. 

Case  1b:  1  Ship,  1  Target,  1  Sector 

Now  consider  the  same  problem  as  Case  la  but  with  a  different  threat  sector, 
as  shown  in  Figure  2.8. 

If  we  divide  the  sector  (a,b)  into  two  sectors  (a,c)  and  (c,b),  we  notice  that 
Eq'n  (2.7)  is  not  valid  in  (c,b).  This  is  so,  because  for  any  target  in  this 
sector,  d*  ,  the  intercept  distance,  is  such  that  the  ship  will  intercept  the 
target  after  it  has  passed  through  the  origin  -  i.e.  after  it  has  passed  the 
carrier.  This  fact  is  illustrated  in  Figure  2.9. 

This  does  not  make  sense  -  there  is  no  reason  to  kill  targets  after  after  they 
have  passed  the  carrier.  In  this  situation,  the  best  that  the  ship  can  do  is 
to  kill  the  target  as  it  is  passing  through  the  origin.  In  other  words,  to 
kill  the  target  at  a  distance  d*-  such  that: 
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2)  Carrier  at  the  Origin  0 

3)  Threat  Sector  (a  ,b  ) ,  a  =  0 

Figure  2.8  Case  1b 
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Legend 

1)  Ship  with  Coordinates  (x‘,y(-) 

2)  Carrier  at  the  Origin  9 

3)  Threat  Sector  (a  ,b  )  Divided  into  Two  Sectors,  (a 

Figure  2.9  Case  1b:  Divided  Sector 


,c  )  and  (c  ,b  ) 
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d*.  =  r. 

> _ • _ 

where  r'  =VJ  x-1  +  y;  is  the  distance  of  ship  i  to  the  carrier. 

Figure  2.10  illustrates  the  result  of  d<y  =  r^. 

In  order  to  include  the  possibility  of  a  Case  la  or  a  Case  lb  sector  we  gener¬ 
alize  d«-  . 

I 

DEFINITION  2.6: 

For  any  ship  i  with  coordinates  (xj  ,y^  ) ,  f*  =  tan  (y-/x)  is  the  angle 
of  ship  i  with  respect  to  the  origin. 

DEFINITION  2.7:  f% 

Let  8  be  the  angle  of  approach  of  a  target.  Then  8  defines  iwo  angles, 
fy  and  .  They  are: 

=  $  -  7/2  and, 
fx  =  8  +  7/2. 

Examples  of  ft  and  are  illustrated  in  Figure  2.11. 

For  any  angle  8  and  any  ship  i  with  coordinates  (x(-  ,y-  ),  d*.  is  the  shortest 
and  therefore  the  perpendicular  distance  from  ship  i  to  the  target.  In  other 
words,  d*.  =  xjsin#  -  yj  cos#  if  and  only  if  f-  $  d.  This  corresponds  to 
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Legend 

1)  8  -  Angle  of  Target  =  60 

2)  ft  »  Angle  of  Ship  i  =  20 

3)  ff  »  8  -  */2  =  -30 

4)  8  +  >r/2  =  150 

Figure  2.11a  Example  of  «,  */?</  ^ 
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saying  that  the  ship  must  he  within  the  shaded  regions  of  the  examples  of  Fig¬ 
ure  2.11.  This  is  expressed  more  formally  in  Definition  2.8. 


DEFINITION  2.8:  GENERALIZED  d* 

i 


d«-  is  the  INTERCEPT  DISTANCE  if  and  only  if 


d«:  =  x-sin0  -  y-cos0 

i  •  i 


■  r  ■ 


otherwise 


For  the  scenario  in  Figure  2.10  we  conclude  that: 

■C 


Prob(ship  i  kills  the  AAW  target) 


f 

*  \  f(0)  x  P_  x 

Ja  A; 


'a 

.b 


-(x;  sine-y^  cosQ) 

e  ^A«  de 

(2.8) 


Xf,9> *  p<*. 


-r 


x  e  dO 

In  this  report,  we  wish  to  examine  only  scenarios  resulting  from  the  type  of 
sector  found  in  Case  la.  This  leads  to  the  followirfg  definition. 
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DEFINITION  2.9:  PERMISSIBLE  SECTOR 


For  a  given  ship  i  with  coordinates  (x* ,  y.  ) ,  a  sector  defined  by 
(a,  b)  is  called  a  PERMISSIBLE  sector  if  and  only  if 


d«;  *  x-  sin0  -  y-  cos0 

»  i  i 


for  all  9 ,  a  v<  8  $  b 


COMMENT  ON  PERMISSIBLE  SECTORS: 


Whether  or  not  a  sector  is  permissible  depends  on  both  the  position  of  the 
ship(s)  as  well  as  the  nature  of  the  sector  (i.e.  size  and  location).  An 
example  of  a  permissible  sector  is  given  in  Figure  2.12a.  The  same  sector  with 
ship  i  in  another  position  is  not  permissible  and  is  shown  in  Figure  2.12b. 
From  now  on  we  study  only  permissible  sectors. 


Case  2:  M  Ships,  1  Target,  1  Threat  Sector 


We  return  to  the  case  of  permissible  sectors  and  consider  Figure  2.13. 

For  the  case  of  M  ships  defending  the  carrier  ,  we  obtain  the  following 
expression: 

Prob(carrier  survives)  =  Prob(at  least  one  ship  kills  the  AAW  target) 


=  1  -  Prob(no  ship  kills  the  AAW  target) 


f  (8)  d  -  nu-p.-  ))  de 


(2.9) 
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Legend 

1)  Threat  Sector  (0  ,120  ) 

2)  6  -  Angle  of  Target  =  120 

3)  f.*  Angle  of  Ship  i  =  60 

4)  fx  =  9  -  ir/2  =  30 

5)  fx  =  9  +  ff/2  =  210 


Target,  1  Sector 

Problem  Statement  and  Mathematical  Formulation 


PAGE  50 


Thus,  we  assume  that  each  of  the  M  ships  will  attempt  to  intercept  the  target 
independently. 


Case  3:  M  Ships,  N  Independent  and  Identically  Distributed  Targets, 
1  Threat  Sector 


Targets  approach  the  carrier  independent  of  each  other.  It  is,  therefore,  easy 
to  generalize  this  problem  for  a  situation  where  the  carrier  is  threatened  by 


more  than  one  target.  We  have  then,  for  the  case  of  N  targets: 


Prob (carrier  survives)  ■ 


-H 


(l  -  Il(l-P  ))  d© 

A, 


(2.10) 


Figure  2.14  shows  the  situation  in  Case  3. 


Case  4:  M  Ships,  N  Independent  and  Identically  Distributed  AAW 
Targets,  0  Independent  and  lndenti cal  1 y  Di stri buted  ASM 
Targets,  2  Identical  Threat  Sectors  (1  A AW;  1  ASM) 


Finally,  we  include  the  possibility  of  both  Kinds  of  threats  -  AAW  and  ASW. 


Figure  2.15  pictures  the  situation  in  Case  4. 
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Prob(carrier  survives)  = 


x 


(2.11) 


2.3  CONCLUDING  REMARKS 


In  this  chapter,  we  have  defined  all  the  necessary  carrier  survival  probabili 


1)  Ship  i  with  Coordinates  (x’.y-) 

i  J  j  ' 

2)  Carrier  at  the  Origin  £ 

3)  Two  Identical  Threat  Sectors:  AAW  in  (a  ,b  )  and  ASW  in  (a  ,b 


Figure  2.15  Case  4:  M  Ships,  N  AM  Targets ,  0  ASW  Targets,  Two  Identical 

Threat  Sectors  (One  AAW,  One  ASW) 
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CHAPTER  3  SINGLE  SHIP  -  SINGLE  SECTOR  AND  DUAL  CAPABILITIES  CASE  STUDIES 


3.0  SUMMARY 


The  purpose  of  this  chapter  is  to  explain,  with  numerical  experimental  evi¬ 
dence,  the  properties  of  the  mathematical  model  developed  in  Chapter  2.  The 
experiments  are  restricted  to  a  scenario  in  which  only  one  ship  is  defending 
the  carrier.  In  the  first  part  of  this  chapter,  the  ship  is  cnaracterized  by 
the  fact  that  it  has  only  AAW  (or  ASW)  capabilities.  In  the  remaining  part  of 
the  chapter,  the  problem  and  the  experiments  are  expanded  to  include  a  ship 
with  dual  (i.e.  AAW  and  ASW)  capabilities,  so  that  it  can  defend  the  carrier 
against  both  types  of  threats. 


We  study  the  effects  on  the  survival  probability  of  the  carrier  and  on  the 


optimum  location  of  the  ship  when  the  size  of  the  threat  sector  is  changed  and 
when  we  introduce  the  possibility  of  multiple  targets.  This  is  all  accom¬ 
plished  through  the  use  of  a  numerical  optimization  program.  Unfortunately, 


even  under  the  simplest  circumstances,  the  optimal  ship  location  problem  can¬ 


not  be  solved  analytically  by  closed-form  formulae. 
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We  explore  the  existence  of  local  maxima.  We  draw  conclusions  about  the 
behaviour  of  the  model  and  set  the  stage  for  the  experiments  of  Chapter  4  - 
the  multiple  ship  case. 


3.1  PROBLEM  STATEMENT  -  SINGLE  SECTOR,  SINGLE  WARFARE  CAPABILITY 


The  purpose  of  the  problem  is  to  maximize  the  probability  of  survival  of  a 
carrier  given  a  single  air  or  sub  threat.  The  threat  is  uniformly  distributed 
over  a  Known  sector.  A  single  AAW  or  ASW  ship  is  defending  the  carrier.  The 
solution  to  the  problem  is  the  optimum  location  of  the  AAW  or  ASW  ship  with 
respect  to  the  carrier;  that  location  which  maximizes  the  survival  probability 
of  the  carrier.  We  formulate  the  problem  in  terms  of  an  AAW  ship;  the  formu¬ 
lation  for  an  ASW  ship  is  analogous. 

The  purpose  of  this  chapter  is  to  perform  parametric  analyses;  to  study  the 
quantitative  impact  of  different  parameter  values  in  our  model  on  the  survival 
probability  of  the  carrier  and  on  the  optimum  location  of  the  ship. 


3.2  MATHEMATICAL  FORMULATION  -  SINGLE  SECTOR,  SINGLE  CAPABILITY 


Hack 
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We  recall  that  we  defined  the  optimum  interception  distance,  the  intercept 
distance,  to  be  the  shortest  distance  between  ship  i  and  the  target 
trajectory.  It  is  at  this  distance  that  the  ship  will  try  to  kill  the  target. 


We  showed  in  Section  2.2.1  that  d»y  was  the  intercept  distance  and  that 


d*.  =  xj  sin#  -  yj-cos0 


(3.1) 


where , 


(x^y,*)  are  the  coordinates  of  ship  i,  and 
8  is  the  angle  of  the  approaching  AAW  target. 


3.2.2.  THE  EFFECT  OF  THE  SHIP'S  DISTANCE  TO  THE  CARRIER 


In  Section  2. 2. 2.1  we  established  that  the  ships  needed  to  maintain  a  certain 
distance  from  the  carrier  in  order  to  detect  the  targets.  We  suggested  incor¬ 
porating  this  constraint  into  the  model  by  introducing  the  function  P^.  •  for 
every  ship  i.  We  recall  from  Eq'n  (2.4)  that 
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p0*  =  1  ■  e 


”(ltA«  r,'  > 


(3.2) 


where , 


rj  =Nxj  +  y(-  is  the  distance  of  ship  i  to  the  carrier,  and 

k^  is  a  scale  factor. 


3.2.3  PROBABILITY  OF  INTERCEPTION 


We  recall  that  each  ship  i  has  associated  with  iu  a  function  P^(-  that  repres¬ 
ents  its  ability  to  kill  AAW  targets.  P^.  was  described  fully  in  Section 
2. 2. 2. 2  ,  Eq'n  (2.5)  and  is  given  by: 


pa;  =  po\,  e 


(-d».  /xA. ) 


where. 


-(kA-  r  ) 

PfA;  =  i  -  e 

d*'  =  (xj sin0  -  jtcostf), 

I is  the  AAW  range  factor  and, 
d  is  the  angle  of  approach  of  the  AAW  target. 
Figure  3.2  shows  P^. 


(3.3) 


(3.2) 

(3.1) 
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3.2.4  PROBABILITY  DISTRIBUTION  OF  THE  TARGET 


The  targets 
2. 2. 2. 3,  let 
Then  f id)  is 


are  uniformly  distributed  over  a  known  sector.  As  in  Section 
fid)  be  the  probability  density  function  of  a  single  target, 
given  by 


fid)  -  l/(b  -  a) 
=  0 


a  ^  d  4  b 
otherwise 


3.2.5  SURVIVAL  PROBABILITY  OF  THE  CARRIER 


We  showed  in  Section  2. 2. 4.1  Case  1  that  the  survival  probability  of  the  car¬ 
rier  was  given  by 


Prob(carrier  survives)  =  Prob(  AAW  ship  kills  the  target  over  sector  (a,b)) 


Figure  3.3  shows  the  area  of  interest  in  the  survival  probability  of  the  car¬ 
rier. 
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3.3  AN  AAW  SHIP  -  EXPERIMENTAL  RESULTS 

The  model  is  complete.  We  will  do  a  sensitivity  analysis  in  order  to  better 
understand  the  properties  of  the  model  and  its  physical  limitations.  We  will 
start  by  investigating  the  effect  of  different  of  and  I  on  the  objec¬ 

tive  function  and  on  the  final  position  of  the  ship.  We  will  examine  the 
results  of  changing  the  size  of  the  sector,  study  the  effect  of  multiple  tar¬ 
gets,  and  explore  the  possibility  of  local  maxima. 

We  will  use  the  following  notation  throughout  the  entire  chapter. 

NOTATION 

The  maximum  probability  of  the  carrier  survival  = 

This  is  the  optimal  value  of  our  objective  function. 

The  optimal  distance  of  the  ship  from  the  carrier -is  denoted  R,  and  its  opti¬ 
mal  angle  is  Thus,  (R  ,}  )  denotes  the  optimal  ship  location  in  polar  coor¬ 
dinates  with  respect  to  the  carrier. 
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3.3.1  SUITABLE  Ka.  AND  £a,  VALUES 


We  consider  again  Eq'n  (3.3)  which  represents  the  AAW  ship's  ability  to  kill 


:argets. 


P4'.  =  e 


(_d\  /^A.  } 


where. 


■ 1  - e 


-<**.•  r,  1 


d«.  =  xj  sin0  -  y ■  cos# 

Z^.  is  the  AAW  range  factor,  and 
6  is  the  angle  of  approach  of  the  target. 


(3.3) 


(3.2) 

(3.1) 


We  want  to  maximize  the  probability  of  survival  of  the  carrier  given  an  AAW 
threat. 


The  optimization  of  this  survival  probability  depends  on  the  ability  of  a  par¬ 
ticular  AAW  ship  to  kill  the  target.  Each  AAW  ship  has  associated  with  it  a 
characteristic  Z^;  value  and  k^.  value.  In  this  section,  we  state  criteria 
for  optimality  (i.e.  >  .9  and  R  =  20)  and  look  for  a  "standard"  AAW  ship: 
one  with  Z^.  and  k^-  which  satisfy  these  conditions  for  optimality. 


What  are  suitable  values  for  k^  and  Z^.  ?  In  order  to  determine  them,  we  ran 
some  experiments  on  the  computer. 
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EXPERIMENT  3.1 


Purpose:  To  find  values  for  and  Z^.  such  that  the  final 

probability  of  survival  of  the  carrier,  Pp,  is  greater  than  .9 
and  R,  the  final  position  of  the  AAW  ship  is  about  20. 

The  sector  extends  from  0-45  and  has  a  single  target  threat. 

Table  3.10  shows  the  value  of  the  fixed  parameters  in  Experiment  3.1. 
Table  3.11  shows  the  effect  of  setting  =  0.2  and  varying  1^.  . 

Table  3.12  shows  the  effect  of  setting  Z^j  *  1000  and  varying  . 

DISCUSSION  OF  RESULTS  ( Table  3.11) 

Consider  Eq’n  (3.4). 

Prob(carrier  survives)  =  \  fit))  x  PA  d#  (3.4) 

Jo 

where. 
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Number 

of 

Ships 

Type 

of 

Ship (s) 

Number 

of 

Sectors 

Type 

of 

Sectors 

Sector 

Number 

of 

Targets 

1 

AAW 

1 

AAW 

0-45 

l 

Table  3.10.  Fixed  Parameters  in  Exp't  3.1:  1  Ship 


1000 

800 

500 

200 

50 

Surv.  Prob.  ^ 

.964 

.959 

.944 

.902 

.794 

Ship  Range  R 

22.5 

21.69 

19.86 

16.43 

11.76 

Ship  Angle  ^ 

22.5 

22.5 

22.5 

22.5 

22.5 

Table  3.11  Effect  of  Varying 

(k^  =0.2) 

. 

xa; 

.05 

0.1 

0.2 

0.3 

0.4 

Surv.  Prob.  P^ 

.815 

.913 

.964 

.979 

.986 

Ship  Range  R 

49.85 

34.49 

22.5 

17.26 

14.14 

Ship  Angle  f 

22.5 

22.5 

22.5 

22.5 

22.5 

Table  3.12  Effect  of  Varying  kA-  ( 1 A.  -  1000) 
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PA-  =  P-  e 

A, 


(-d«  /i  ) 

i  ** « 


and. 


P.  =  1  -  e 

Or, 


-(1ca;  r,  } 


d*.  =  x-  sin#  -  y-  cos# 

i  •  • 

It  is  evident  that  as  increases,  the  probability  of  survival  increases. 

So  we  expect  in  Table  3.11,  that  a  larger  value  would  yield  a  larger  Pp  . 

This  seems  to  be  the  case.  We  would  also  suspect  that  a  larger  1^-  would 

result  in  the  distance  of  the  ship  from  the  carrier,  R,  to  be  larger  because  a 

bigger  range  means  that  a  ship  need  not  be  so  close  to  the  carrier.  Because  of 

the  symmetry  of  the  problem,  and  the  uniform  distribution  of  the  target  within 

o 

the  threat  sector,  the  final  angle  of  the  ship  in  all  the  examples  is  22.5. 


DISCUSSION  OF  RESULTS  (Table  3.12) 

Again,  from  observing  Eq'n  (3.4),  we  see  that  as  the  value  of  )c^(  increases, 
the  value  of  the  probability  of  survival  increases.  This  is  consistent  with  f^. 
increasing  with  k  . 


3. 3. 2. 2  A  CHANGE  IN  SECTOR  SIZE 
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We  now  increase 


size  of  the  sector  from  0-90  and  consider  the  next 


experiment. 

EXPERIMENT  3.2 

Purpose:  To  increase  the  sector  size  from  0-90  and,  using  the  same 

values  of  k^  and  as  in  Experiment  3.1,  compare  the  values 

of  P_p  and  P  with  those  of  Experiment  3.1. 

Table  3.20  shows  the  values  of  the  fixed  parameters  in  Experiment  3.2. 

Table  3.21  shows  the  effect  of  setting  kft-  =  0.2  and  varying  . 

Table  3.22  shows  the  result  of  setting  S^;  =  1000  and  varying  k  . 

DISCUSSION  OF  RESULTS  (Table  3.21  and  Table  3.22' 

As  we  expect,  when  the  threat  sector  is  expanded  the  ship  is  no  longer  as 
effective  and  the  carrier  survival  probability  value,  P^  ,  drops.  The  optimal 
ship-to-carrier  distance,  R,  decreases  because  the  ship  must  be  pulled  toward 

the  carrier  in  order  to  better  cover  the  carrier.'  Once  more,  because  of  sym- 

o 

metry,  the  angle  of  the  defending  ship  is  45,  at  the  center  of  the  threat 
sector. 

In  the  next  experiment,  we  observe  what  happens  when  we  change  the  size  of  the 
sector  and  leave  all  other  parameters  fixed. 
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Number 

of 

Ships 

Type 

of 

Ship (s) 

Number 

of 

Sectors 

Type 

of 

Sector (s) 

Sector  Number 

of 

Targets 

i 

AAW 

1 

AAW 

0-90  1 

Table  3.20  Fixed  Parameters  in  Exp't  3.2:  1  Ship 


1000 

800 

500 

200 

50 

Surv,  Prob.  P^. 

.918 

.906 

.877 

.802 

.644 

Ship  Range  R 

17.57 

16.75 

15.07 

12.02 

8.09 

Ship  Angle  f 

45 

45 

45 

45 

45 

Table  3.21  Effect  of  Varying 

(k..  »  0.2) 

A, 

• 

2  a; 

.05 

0.1 

0.2 

0.3 

0.4 

Surv.  Prob.  P^ 

.672 

.823 

.918 

.951 

.966 

Ship  Range  R 

34.67 

25.47 

17.57 

13.77 

11.47 

Ship  Angle  f 

45 

45 

45 

45 

45 

Table  3.22  Effect  of  Varying  k ^  ( IA.  =»  1000 ) 

Chapter  3  Single  Ship  -  Single  Sector  and  Dual  Capabilities  Case  Studies 


EXPERIMENT  3.3 


Purpose:  To  observe  the  effect  on  the  carrier  survival  probability,  P  , 
and  on  the  ship  range,  R,  when  we  vary  the  size  of  the  sector. 


Table  3.30  shows  the  values  of  the  fixed  parameters  in  Experiment  3.3. 
Table  3.31  shows  the  effect  of  setting  =  1000  and  k^.  =  .2  and  varying 


the  size  of  the  sector. 


DISCUSSION  OF  RESULTS  (Table  3.31) 

Again,  we  observe  that  as  the  size  of  the  sector  is  increased,  it  becomes  more 
difficult  for  the  ship  to  protect  the  carrier  and  as  a  result,  the  optimal 
survival  probability,  Pp  ,  decreases.  Also,  as  the  size  of  of  the  sector 
increases,  the  optimal  ship-to-carrier  distance,  R,  shrinks  because  the  ship 
needs  to  move  closer  to  the  carrier  in  order  to  better  defend  it. 


3. 3. 2. 3  MULTIPLE  TARGETS 


What  happens  to  the  probability  of  survival  and  the  final  location  of  the  ship 
when  the  number  of  targets  in  a  sector  is  increased  to  N? 

The  survival  probability  of  the  carrier  based  on  this  scenario  is: 


Prob(carrier  survives) 


a; 


(3.5) 
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Number  Type  Number  Type  fcA. 

of  of  of  of 

Ships  Ship(s)  Sectors  Sectors 

1  AAW  1  AAW  1000  0.2 


Table  3.30  Fixed  Parameters  in  Exp't  3.3:  1  Ship 


Sector 

0-15 

0-20 

0-45 

0-60 

0-75 

0-90 

0-105 

0-120 

Pf 

.992 

.988 

.964 

.948 

.933 

.918 

.904 

.891 

R 

31.63 

29.18 

22.5 

20.40 

18.80 

17.57 

16.59 

15.81 

t 

7.5 

10 

22.5 

30 

37.5 

45 

52.5 

60 

Table  3.31  Effect  of  Varying  the  Sector  Size 
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We  expect  to  decrease  since  the  threat  is  more  dense.  Since  the  threats  are 
independent  and  identically  distributed,  we  expect  that  R,  the  optimal 
ship-to-carrier  distance,  will  not  change. 

CLAIM  3.1: 

For  independently  and  identically  distributed  targets,  R 
does  not  change  with  the  number  of  targets. 

» 

Proof 

Please  refer  to  Appendix  A  for  the  proof. 

We  now  verify  this  with  experimental  results. 

EXPERIMENT  3.4 

Purpose:  To  see  the  effect  of  multiple  targets  on  and  R. 

Table  3.40  shows  the  values  of  the  fixed  parameters  in  Experiment  3.4. 

Table  3.41  shows  the  result  of  varying  the  number  of  targets. 
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Number 

of 

Ships 

Type 

of 

Ship (s) 

Number 

of 

Sectors 

Type  Sector 

of 

Sector 

2a. 

** 

1 

AAW 

1 

AAW  0-45 

1000 

0.2 

Table  3.40  Fixed  Parameters  in  Exp't  3.4:  1  Ship 


Number  of  Targets 

1 

2 

5 

10 

20 

Surv.  Prob.  P^ 

.964 

.930 

.834 

.696 

.485 

Ship  Range  R 

22.58 

22.58 

22.58 

22.58 

22.58 

Ship  Angle  f 

22.5 

22.5 

22.5 

22.5 

22.5 

Table  3.41  Effect  of  Varying  the  Number  of  Targets 
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DISCUSSION  OF  RESULTS  (Table  3.41) 

We  see  that  the  ship-to-target  distance,  R,  remains  constant  (as  does  the 
optimal  angle  of  the  ship)  as  the  number  of  targets  in  the  sector  changes. 
This  is  a  direct  result  of  CLAIM  3.1  above. 

We  also  observe  that  the  carrier  survival  probability,  Pp  ,  decreases  rapidly 
as  the  number  of  targets  increases  substantially.  When  we  double  the  sector 


size 

to  0-90, 

and  set  I =  1000 

a  )C^F  0  •  2  r 

then  =  .918  (Table 

3.21)  . 

In 

Experiment  3.3, 

Table  3.31,  with 

I*  =  1000 

,  k.  =  0.2  and  a  0-45 

A,'  • 

sector 

but 

with 

double  the 

concentration  of 

targets. 

Pp  =  .930.  Comparing 

these 

two 

results,  it  is  more  detrimental  in  this  case  to  double  the  sector  size  than  to 
double  the  concentration  of  targets. 


3.3.3  EXISTENCE  OF  LOCAL  MAXIMA 


In  this  section,  we  explore  the  possibility  of  local  maxima.  That  is,  we  vary 
the  initial  position  of  the  AAW  ship  to  determine  whether  it  has  any  influence 
on  the  optimal  position  of  the  ship.  If  the  optimal  location  of  the  ship  is 
independent  of  the  initial  positions, we  conclude  that  there  are  no  local 
maxima.  We  make  intelligent  choices  for  the  initial  conditions.  That  is,  we 
choose  a  position  which  seems  "reasonable".  If  we  want  the  optimal 
ship-to-carrier  distance,  R,  to  be  about  20,  we  would  not  start  with  an  ini- 
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tial  position  of  the  ship  at  (80,80).  Not  only  may  it  he  costly  because  of  the 
incremental  number  of  iterations  that  might  be  required  to  converge  to  the 
(desired)  optimal  solution,  but  also  such  choices  increase  the  probability  of 
falling  on  a  local  maximum,  which  yields  an  undesirable  (i.e.  too  far  away) 
location.  In  some  cases,  choosing  an  initial  condition  which  is  "far"  from  the 
'~-;»rrier  (such  as  (80,80)  in  Table  3.52)  results  in  the  same  ^  but  shifts  the 
optimal  angle  of  the  ship  by  180  degrees.  These  are  limitations  of  the  model. 


We  recall  that  we  are  dealing  with  permissible  sectors.  A  permissible  sector 
was  defined  (Def’n  2.8)  to  be  one  in  which,  for  a  given  ship  with  coordinates 


(x;  ,yj-  ) ,  d*  ,  the  intercept  distance,  is  given  by: 


d*  =  x‘sin0  -  y.cosfl 

i  <  * 


for  all  d,  a  $  d  f  b 


(See  Section  2. 2. 3.1,  Figure  2.12b  for  an  example  of  a  position  (or  initial 
condition)  of  a  ship  which  does  not  generate  a  permissible  sector.) 


The  initial  condition  that  we  select  for  our  ship  must  generate  a  permissible 
sector. 
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EXPERIMENT  3.5 


Purpose:  To  test  for  the  existence  of  local  maxima.  We  keep  all  other 
parameters  fixed  and  vary  the  initial  position  of  the  ship. 

Table  3.50  shows  the  values  of  the  fixed  parameters  in  Experiment  3.5. 

Table  3.51  shows  the  result  of  changing  the  initial  position  of  the 
ship. 

DISCUSSION  OF  RESULTS  (Table  3.51} 

we  observe  that  for  "reasonable"  initial  conditions  (i.e.  ones  that  generate 
permissible  sectors  and  are  not  too  far  from  the  desired  optimal 
ship-to-carrier  distance,  R,  the  model  appears  not  to  be  sensitive  to  changes 
in  the  initial  condition.  In  our  examples,  at  worst  the  optimal  angle  of  the 
ship  shifted  by  180  degrees  and  this  presented  an  unrealistic  solution  for  our 
problem.  The  optimal  survival  probability,  I>,  did  not  change. 


3.4  AW  ASW  SHIP  -  EXPERIMENTAL  RESULTS 


In  this  section,  we  perform  experiments  on  an  ASW  ship  similar  to  those  that 
we  performed  on  the  AAW  ship  in  Section  3.3.  The  purpose  is  to  discover  any 
trends  in  the  optimal  carrier  survival  probability,  P^.  ,  and  in  the  optimal 
position,  R,  of  the  ASW  ship. 
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Number 

Type 

Number 

Type 

Sector 

1a; 

-A; 

Number 

of 

of 

of 

of 

of 

Ships 

1 

Ship(s) 

AAW 

Sectors 

1 

Sector 

AAW 

0-45 

1000 

0.2 

Targets 

1 

Table  3.50  Fixed  Parameters  in  Exp't  3.5:  1  Ship 


Init.  Cond'n 

(5,1) 

(20,20) 

(1,30) 

(80,80) 

(40,0) 

(0,3) 

(1,-1) 

Surv.  Prob. 

.964 

.964 

.964 

.964 

.964 

.964 

.964 

Ship  Range  R 

22.58 

22.58 

22.58 

22.58 

22.58 

22.58 

22.58 

Ship  Angle  ^ 

22.5 

22.5 

22.5 

202.5 

22.5 

22.5 

22.5 

Table  3.51  Effect  of  Changing  the  Initial  Condition 
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We  will  discuss  the  results  of  changing  the  size  of  the  sector,  study  the 
effect  of  multiple  targets,  and  review  the  possibility  of  local  maxima. 

The  reader  should  be  reminded  of  the  fact  that  in  our  model  the  motion  and 
interception  of  submarines  is  treated  in  a  similar  manner  as  that  of  an  air 
threat.  Each  submarine  target  is  assumed  to  follow  a  straight  line  path  to  the. 
carrier.  This  assumption  implies  that  we  are  treating  a  submarine  threat  as 
though  it  has  to  closely  approach  the  carrier  before  it  launches  a  torpedo;  it 
does  not  attempt  to  evade  ASW  weapons.  Also,  the  fact  that  submarines  can 
create  air  threats,  by  launching  missiles  against  the  carrier,  is  not 
addressed  (unless  separately  modeled  in  the  AAW  problem)  in  our  model*  . 

We  emphasize  that  this  "symmetry"  of  AAW  and  ASW  threats  is  a  major  limitation 
in  the  practical  application  of  our  model.  Further  research  is  needed  to 
improve  the  relevance  of  the  ASW  part  of  our  model. 


3.4.1  SUITABLE  K. 


AND  2^; 


VALUES 
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We  consider  an  analogous  equation  to  Eq'n  (3.3)  -  the  ASW  ship's  ability  to 
kill  ASW  targets. 


<V  /zs  > 


where, 


) 


(3.5) 


d»  =  x;sin$  -  yrcostf 

i  *  1 

X5;  is  the  ASW  range  factor,  and 
d  is  the  angle  of  approach  of  the  ASW  target. 


Our  purpose  is  to  maximize  the  probability  of  survival  of  the  carrier,  given  a 
single  ASW  target  threat.  The  optimization  of  this  survival  probability 
depends  on  the  ability  of  a  particular  ASW  ship  to  kill  the  target.  Each  ASW 
ship  has  associated  with  it  a  characteristic  Xs-  and  k^.  value.  We  impose 
criteria  for  optimality:  the  value  of  must  exceed  .9  and  the  optimal 

ship-to-carrier  distance  R,  must  be  approximately  SO. 

We  look  for  an  ASW  ship  that  satisfies  these  conditions,  i.e.  we  want  to 
determine  the  characteristic  k^  and  X^-  values  of  an  ASW  ship. 

What  are  suitable  values  for  k^  and  X^.  ?  We  run  the  next  experiment  to 
determine  them. 
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EXPERIMENT  3.6 

Purpose:  To  find  values  for  k5-  and  1^-  such  that  the  final 

probability  of  survival  of  the  carrier,  ,  is  greater  than  .9 
and  R,  the  final  position  of  the  AAW  ship  is  about  50. 

The  sector  extends  from  0-45  and  has  a  single  target  threat. 


Table  3.60  shows  the  value  of  the  fixed  parameters  in  Experiment  3.6. 
Table  3.61  shows  the  effect  of  setting  k«..  =  0.07  and  varying  Z<  . 

J| 

Table  3.62  shows  the  effect  of  setting  Z5;  =  2000  and  varying  . 
Table  3.63  shows  the  effect  of  setting  =  8000  and  varying  k^.  . 


DISCUSSION  OF  RESULTS  (Table  3.61) 

Consider  Eq'n  (3.6). 


Prob (carrier  survives)  = 


x 


where. 


(3.6) 
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Number 

Type 

Number 

Type 

Sector 

Number 

of 

of 

of 

of 

of 

Ships 

Ship 

Sectors 

Sector 

Targets 

1 

ASW 

1 

ASW 

0-45 

1 

Table 

3.60  Fixed  Parameters  in 

Exp't  3.6: 

1  Ship 

8000 

5000 

3000 

2600 

2200 

2000 

1000 

Pf 

.964 

.950 

.931 

.925 

.917 

.912 

.869 

R 

64.30 

58.99 

53.38 

51.84 

50.06 

49.05 

42.00 

t 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

Table  3.61  Effect  of  Varying  1^.  (k^.  =  .07) 
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ra 


V-;' 
►“ . 


** 

.02 

in 

o 

• 

.07 

.08 

.09 

.10 

p-f 

.692 

.871 

.912 

.925 

.935 

.943 

R 

91.14 

59.04 

49.05 

45.40 

42.33 

39.72 

t 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

Table  3.62  Effect  of  Varying  k<-  ^1, 


2000) 


H- 

.02 

.05 

.07 

00 

o 

• 

.09 

.10 

P* 

.837 

.943 

.964 

.969 

.974 

.978 

R 

132.63 

79.43 

64.30 

58.95 

54.53 

50.81 

t 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

Table  3.63  Effect  of  Varying  (ls.  -  8000) 
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ps;  *  ?<-s;  - 


(-a.,  /i  ) 

I  J) 


Pr  *  1  -  e 
*  ^  ■ 


-(it.,  r- 

V|  I 


=  Xj  sin# 


y  cos# 
» 


It  is  evident  that  as  increases,  the  probability  of  survival  increases. 

So  we  expect  in  Table  3.61,  that  a  larger  1^-  value  would  yield  a  larger  Pp  . 
This  seems  to  be  true.  We  would  also  expect  that  a  larger  would  result  in 

the  distance  of  the  ship  from  the  carrier,  R,  to  be  larger  because  a  bigger 
range  means  that  a  ship  need  not  be  so  close  to  the  carrier.  Because  of  the 

symmetry  of  the  problem,  and  the  uniform  distribution  of  the  target,  the  final 

o 

angle  of  the  ship  in  all  the  examples  is  22.5. 


DISCUSSION  OF  RESULTS  (Tables  3.62  and  3.63) 

Again,  from  observing  Eq'n  (3.6),  we  see  that  as  the  value  of  increases, 
the  value  of  the  probability  of  survival  increases.  This  is  consistent  with  Pp 
increasing  with  X^-  .  A  larger  increases  the  Pp  value  and  also  increases 
the  value  of  R. 
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3. 3. 2. 2  A  CHANGE  IN  SECTOR  SIZE 


We  now  increase  the  size  of  the  sector  from  0-90  and  consider  the  next 


experiment. 


EXPERIMENT  3.7 


Purpose:  To  increase  the  sector  size  from  0-90  and,  using  the  same 

values  of  kg.  and  Ig.  as  in  Experiment  3.6,  compare  the  values 
of  Pr  and  R  with  those  of  Experiment  3.6. 


Table  3.70  shows  the  values  of  the  fixed  parameters  in  Experiment  3.7. 

Table  3.71  shows  the  effect  of  setting  k<-.  =  0.07  and  varying  I,.  . 

Table  3.72  shows  the  result  of  setting  Zgj  =  2000  and  varying  kg. 

DISCUSSION  OF  RESULTS  (Table  3.71  and  Table  3.72 1 

As  we  expect,  when  the  threat  sector  size  is  increased  the  ship  can  no  longer 
defend  the  carrier  as  well,  and  the  value  of  P^  decreases.  Also,  the  optimal 
ship-to-carrier  distance,  R,  decreases  since  the  ship  must  get  closer  to  the 
carrier  in  order  to  better  defend  it. 
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Number 

Type 

Number 

Type  Sector 

Number 

of 

Ships 

of 

Ship 

of 

Sectors 

of 

Sector 

of 

Targets 

1 

ASW 

1 

ASW  0-90 

l 

Table  3.70  Fixed  Parameters  in  Exp't  3.7: 

1  Ship 

** 

8000 

5000 

3000 

2600 

2200 

2000 

1000 

pf 

.917 

.890 

.855 

.843 

.829 

.821 

.752 

R 

49.99 

45.11 

40.03 

38.66 

37.08 

36.20 

30.09 

t 

45 

45 

45 

45 

45 

45 

45 

Taole  3.71  Effect  of  Varying  X* 3|-  ( k -  .07) 


ks; 

.02 

in 

o 

• 

.07 

.08 

.09 

.10 

pf 

.522 

.754 

.821 

.844 

.862 

.877 

R 

60.22 

42.30 

36.20 

33.87 

31.87 

30.14 

t 

45 

45 

45 

45 

45 

45 

Table  3.72  Effect  of  Varying  k (Ic.  -  2000) 

3 1 
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3. 4. 2. 3  MULTIPLE  TARGETS 


What  happens  to  the  probability  of  survival  and  the  final  location  of  the  ship 
when  the  number  of  targets  in  a  sector  is  increased  to  Q  ? 

The  survival  probability  of  the  carrier  is: 


Prob(carrier  survives) 


s; 


(3.5) 


The  situation  is  completely  analogous  to  the  one  in  Section  3. 3. 2. 3  for  AAW 
ships,  and  we  expect  the  same  results.  We  expect  P^  to  decrease  since  the 
threat  is  more  dense.  Since  the  threats  are  independent  and  identically  dis¬ 
tributed,  we  expect  that  R,  the  optimal  ship-to-carrier  distance,  will  not 
change . 

We  present  the  next  experiment  to  support  this  hypothesis. 
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EXPERIMENT  3.8 

Purpose:  To  see  the  effect  of  multiple  targets  on  Pp  and  R. 

Table  3.80  shows  the  values  of  the  fixed  parameters  in  Experiment  3.8. 
Table  3.81  shows  the  result  of  varying  the  number  of  targets. 


DISCUSSION  OF  RESULTS  (Table  3.811 

For  a  maximum  of  twelve  targets,  the  program  performs  well.  The  value  of  ^ 
dropped  as  the  number  of  targets  increased,  as  we  expected.  The  position  of 
the  ship  was  not  expected  to  change,  (based  on  CLAIM  3.1  of  Section  3. 3.2. 3), 
and  it  did  not.  Why  can  the  program  not  worlc  for  a  number  of  targets  greater 
than  fifteen?  There  is  no  relative  function  convergence,  only  convergence  of 
the  iterates,  (xj ,y- ) ,  and  this  occurs  in  two  iterations.  This  has  no  physical 
importance.  Because  we  are  solving  this  problem  numerically,  using  a  standard 
optimization  routine,  we  are  at  the  mercy  of  its  limitations.  The  result  is 
that  we  must  put  an  upper  limit  oi  fifteen  on  the  number  of  targets.  However, 
such  limitations  occur  only  when  the  carrier  survival  probability  is  almost 
zero;  hence,  such  shortcomings  have  little  impact  on  more  realistic  scenarios. 


The  following  experiment  shows  what  happens  with  different  values  of  )c^.  and 

Z5.' 
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Number 

of 

Ships 

Type 

of 

Ship 

Number 

of 

Sectors 

Tffie 

of 

Sector 

Sector 

U 

*Si 

1 

ASW 

1 

ASW 

0-45 

2000 

.07 

Table  3.80  Fixed  Parameters  in  Exp't  3.8:  1  Ship 


Number 


of 

1 

2 

5 

10 

12 

14 

15 

20 

Targets 

-10 

pf 

.912 

.832 

.633 

.401 

.334 

.3.01 

.1  X  10 

.4  X  10 

R 

49.05 

49.05 

49.05 

49.05 

49.05 

49.05 

5.09 

5.09 

t 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

11.31 

11.31 

Table  3.81  Effect  of  Varying  the  Number  of  Targets 
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EXPERIMENT  3.9 

Purpose:  To  see  the  effect  of  multiple  targets  on  and  R  when  I  =8000 

and  k..  =.l  . 

J\ 

Table  3.90  shows  the  values  of  the  fixed  parameters  in  Experiment  3.9. 

Table  3.91  shows  the  result  of  varying  the  number  of  targets. 

DISCUSSION  OF  RESULTS  (Table  3.91) 

As  the  number  of  targets  increased  from  one  to  eighteen,  the  value  of 
decreased,  and  the  position  of  the  ship  did  not  change.  This  we  expected  since 
the  ship  is  no  longer  as  effective,  and  the  targets  are  independent  of  each 
other.  When  the  number  of  targets  was  set  to  nineteen,  the  program  failed.  As 
in  Experiment  3.8,  there  was  no  relative  function  convergence,  only  conver¬ 
gence  of  the  iterates,  which  is  meaningless.  It  is  worthwhile  to  note  that 
for  Zg.  =  2000,  kg(.  =  .07,  (Exp't  3.8)  the  program  could  sustain  only  fourteen 
targets.  In  Experiment  3.9,  with  !«.-  =  8000,  kfi.  =  .1,  the  program  worked  for 
nineteen  targets.  Also,  in  Experiment  3.4  of  the  AAW  ship,  the  program  worked 
for  twenty  targets. 


3.4.3  EXISTENCE  OF  LOCAL  MAXIMA 
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In  Section  3.3.3,  we  considered  the  existence  of  local  maxima.  That  is,  we 
observed  the  effect  of  the  initial  position  of  the  AAW  ship  on  its  optimal 
position.  In  this  section,  we  do  the  same.  Again  we  make  intelligent  choices 
for  the  initial  conditions  for  the  same  reasons  presented  in  Section  3.3.3. 


Number 

of 

Ships 

Type 

of 

Ship 

Number 

of 

Sectors 

Type 

of 

Sector 

Sector 

2*; 

1 

ASW 

1 

ASW 

0-45 

8000 

Table  3.90  Fixed  Parameters  in  Exp't  3.9:  1  Ship 


Number 


of 

1 

2 

5 

10 

15 

18 

19 

20 

Targets 

.978 

.956 

.894 

.800 

.716 

.670 

.3  x  10 

.1  x  10 

R 

50.81 

50.81 

50.81 

50.81 

50.81 

50.81 

5.09 

5.09 

t 

22.5 

22.5 

22.5 

22.5 

22.5 

22.5 

11.3 

11.3 

Table  3.91  Effect  of  Varying  the  Number  of  Targets 
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EXPERIMENT  3.10 

Purpose:  To  determine  if  there  are  any  local  maxima. 

Table  3.100  shows  the  values  of  the  fixed  parameters  in  Experiment  3.10. 

Table  3.101  shows  the  effect  on  and  R  of  changing  the  initial 
condition. 

DISCUSSION  OF  RESULTS  (Table  3.101) 

We  observe  that  for  "reasonable"  initial  conditions,  the  model  appears  to  be 
robust  in  the  sense  that  the  optimal  survival  probability,  Pp ,  and  the  optimal 
ship-to-carrier  distance,  R  ,  are  unaffected. 


3.5  PROBLEM  STATEMENT  -  SINGLE  SHIP  WITH  DUAL  CAPABILITIES 


The  purpose  of  the  problem  is  to  maximize  the  probability  of  survival  of  a 
carrier  given  a  single  air  and  a  single  sub  threat.  The  carrier  is  being 
defended  by  either  a  single  AAW  ship  or  a  single  ASW  ship.  The  aaw  and  ASW 
threat  sectors  are  identical.  The  solution  to  the  problem  is  the  optimum 
location  of  the  AAW  (ASW)  ship  with  respect  to  the  carrier?  that  position 
which  maximizes  the  survival  probability  of  the  carrier.  The  problem  is  formu¬ 
lated  in  terms  of  an  AAW  ship;  the  formulation  for  an  ASW  ship  is  analogous. 
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Number 

of 

Ships 

Type 

of 

Shi£ 

Number 

of 

Sectors 

Type 

of 

Sector 

Sector 

Z*; 

Number 

of 

Targets 

1 

ASW 

1 

ASW 

in 

l 

o 

8000 

.10 

l 

Table  3. 

100  Fixed  Parameters  in 

Exp't  3.10: 

1  Ship 

I nit.  Cond'n 

(0,3) 

(30,0) 

(80,80) 

(1,-1) 

surv.  Prob.  Pp 

.978 

.978 

.978 

.978 

Ship  Range  R 

50.81 

50.81 

50.81 

50.81 

Ship  Angle  f 

22.5 

22.5 

22.5 

22.5 

Table  3.101  Effect  of  Changing  the  Initial  Condition 
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3.6  MATHEMATICAL  FORMULATION  -  DUAL  CAPABILITIES 

The  carrier  is  located  at  the  origin.  The  AAW  ship  has  coordinates  (x; ,y;  ) 
which  represent  its  position.  There  are  two  identical  threat  sectors,  (a,b) , 
one  AAW,  one  ASW.  The  scenario  is  pictured  in  Figure  3.4. 

3.6.1  OPTIMUM  INTERCEPTION  DISTANCE 

The  intercept  distance  is  the  shortest  distance  between  ship  i  and  the  target 
trajectory.  In  this  case,  it  is  the  shortest  distance  between  ship  i  and 
either  an  AAW  or  ASW  target.  From  Section  2.2.1,  d*  is  the  intercept  distance 
and 

d>«  =  x^sin#  -  y.cos#  for  all  9,  a  <<  9  b. 

3.6.2  THE  EFFECT  OF  THE  SHIP'S  DISTANCE  TO  THE  CARRIER 

As  in  Section  2. 2. 2.1,  we  have  a  function  that  reflects  the  ships'  need  to 
maintain  a  certain  distance  from  the  carrier.  We  recall  from  Eq'n  (3.4)  that 
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■(H 


PfAi  =  1  "  e 


a; 


(3.4) 


where, 

J  ~  *- 

r-  = lx'  +  y-  is  the  distance  of  ship  i  to  the  carrier,  and 
is  a  parameter. 

In  the  case  of  an  AAW  ship  with  dual  capabilities,  it  is  no  longer  sufficient 
to  represent  the  effect  of  the  ship-to-carrier  distance  by  F^.  of  Sq'n  (3.4) 
This  is  so  because  an  additional  threat  (the  ASW  one)  of  a  different  nature 
has  been  introduced.  The  optimal  position  of  the  AAW  ship  will  depend  on  both 
its  AAW  and  ASW  defense  capabilities.  The  latter  are  captured  independently 
in  the  optimal  ship-to-carrier  distance,  R.  Consequently,  we  introduce  an 
analogous  function  to  that  of  Eq'n  (3.4),  and  say  that  the  effect  of  the  AAW 
ship's  distance  to  the  carrier  is  described  completely  by  both  of  these  func¬ 
tions. 

This  new  function  is: 
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v 


"(ks;  r.  } 

P,  -  1  -  e  ,  k;  >  0  (3.6) 

6i 


where, 

JT - *■ 

ri  =  *xi  +  y,*  is  the  distance  of  ship  i  to  the  carrier,  and 
kc.  is  a  parameter. 

Ji 


3.6.3  PROBABILITY  OF  INTERCEPTION 


The  ability  of  a  ship  to  kill  AAW  targets  is  given  by 

(3.3) 

(3.2) 

(3.1) 

Similarly,  we  said  that  the  ability  of  a  ship  to  kill  ASW  targets  was  given  by 
an  analogous  function,  P^.  ,  where 


pa;  =  p0»; 


where , 


V 1  - e 


(~ d*-  /I%  ) 


-(k.-  r  -  ) 

i 


d*  =  xjsinfl  -  y; cosfl  , 

£^-  is  the  AAW  range  factor  and 
6  is  the  angle  of  approach  of  the  AAW  target. 


4 
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where , 


(-d*X  /2S.  ) 
1 


V 


<3*  =  xjsin#  -  yjcosfl, 

1^.  the  ASW  range  factor  and, 

8  is  the  angle  of  approach  of  the  ASW  target. 


(3.7) 


3.6.4  PROBABILITY  DISTRIBUTIONS  OF  THE  TARGETS 

Both  the  AAW  and  the  ASW  targets  are  uniformly  distributed  over  the  same 
(known)  sector.  Let  f(0)  be  the  probability  density  function  of  a  single  AAW 
target.  Then  f (0)  is  given  by 

f (8)  =  l/(b  -  a)  a  <  8  «  b 

*  0  otherwise 

If  g (8)  is  the  probability  density  function  of  a  single  ASW  target,  then 
g (8)  »  f (8),  since  they  are  distributed  identically. 


3.6.5  SURVIVAL  PROBABILITY  OF  THE  CARRIER 
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The  probability  that  the  carrier  survives  is  the  probability  that  the  AAW  ship 
kills  the  AAW  target  and  the  ASW  target.  This  is  given  by 

Prob (carrier  survives)  *  Prob(AAW  ship  kills  the  AAW  and  the  ASW 

target  in  sector  (a,b)) 

*  1  f(8)  P*.  ad  •  f  f(6)  p-  ae  (3.8) 

Jo  *  Ja  *' 

3.7  AW  AAW  SHIP  -  EXPERIMENTAL  RESULTS  (DUAL  CAPABILITIES) 

We  have  built  a  model  that  captures  the  dual  capabilities  of  an  AAW  ship,  and 
the  survival  probability  of  the  carrier  under  an  AAW  and  an  ASW  threat.  Our 
goal  is  to  maximize  the  probability  of  survival  of  the  carrier  given  an  AAW 
and  an  ASW  threat.  The  optimization  of  this  survival  probability  depends  on 
the  AAW  ship's  ability  to  kill  the  AAW  and  the  ASW  target.  We  impose  criteria 
for  optimality:  the  value  of  the  carrier  survival  probability,  P^  ,  must 
exceed  .80,  and  the  the  optimal  ship-to-target  distance, R,  must  be  approxi¬ 
mately  20  miles.  In  Experiment  3.1  we  imposed  similar  conditions:  we  said 
that  Pf  had  to  exceed  .9  and  that  R  should  be  about  20  miles,  in  that  exper¬ 
iment,  the  AAW  ship  had  only  AAW  capabilities  and  was  defending  the  carrier 
against  a  single  AAW  target.  In  this  section,  we  consider  the  same  ship  but 
with  both  AAW  and  ASW  capabilities,  required  to  defend  the  carrier  against 
both  an  AAW  and  an  ASW  target.  We  do  not  expect  the  ship  to  perform  as  well 
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because  of  the  dual  threat.  We  require  it,  however,  to  remain  about  20  miles 
from  the  carrier. 


We  wish  to  mathematically  describe  the  AAW  ship  that  satisfies  these  criteria 

for  optimality  by  finding  the  numerical  values  for  kv  ,  k..  ,  I  and  X<  .  . 

0\ 


3.7.1  SUITABLE  K a;  ,  K 


VALUES 


We  consider  Eq'ns  (3.3)  and  (3.7)  which  represent  the  AAW  ship's  ability  to 
kill  AAW  and  Asw  targets  respectively. 


p«  ;  *  V  • 


'■*1  /z«,-  > 


where , 


(3.3) 


"<*»;  r;  1 


V  1  - e 


d*  *  x'sinfl  -  y*cos0 

i  i  *1 

2 Ai  the  AAW  range  factor  and, 
d  is  the  angle  of  approach  of  the  AAW  target. 


Also, 
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P*  *  Pr  e 

*v, 


(-<3S  /Z-.  ) 

■ 


where , 


-Oc  .  r  -  ) 

a.  ■  1  -  e  , 

Vi'. 


V 


d*  *  Xjsintf  -  y(-cos0 

Ij.  is  the  ASV  range  factor  and 

8  is  the  angle  of  approach  of  the  ASM  target. 


(3.7) 


What  are  suitable  values  for  kA.  ,  k- .  ,1  .  and  I  ? 

5| 


EXPERIMENT  3.11 

Purposes  To  find  values  for  k^j  ,  k^.  ,  1^.  ,  and  Z^  ,  such  that  the 

probability  of  survival  of  the  carrier,  ,  is  greater  than  .8 
and  R,  the  optimal  position  of  the  AAW  ship  is  about  20. 

The  AAW  and  ASW  sectors  extend  from  0-45  and  both  have  a  single 
target  threat. 

Table  3.110  shows  the  value  of  the  fixed  parameters  in  Experiment  3.11. 

Table  3.111  shows  the  effect  of  setting  Zg.  *  500  and  varying  k^.  . 
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Number 

of 

Ships 

Type 

of 

Ship 

Number 

of 

Sectors 

Type 

of 

Sectors 

Sectors 

Number 

of 

Targets 

Zk  £a; 

1 

AAV 

2 

1  AAW 

0-45 

1  AAW 

1000  .2 

1  ASW 

0-45 

1  ASW 

Table  3.110  Fixed  Parameters  in  Exp't  3.11: 

1  Ship,  Dual 

Capabilities 

** 

500 

500 

500 

500 

.02 

.07 

• 

.09 

.1 

Surv.  Prob. 

.448 

.769 

.817 

.834 

Ship  Range  R 

48.13 

31.70 

28.57 

27.34 

Ship  Angle  6 

22.5 

22.5 

22.5 

22.5 

Table  3.111  Effect  of  Varying  (Ii;  -  500/ 
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DISCUSSION  OF  RESULTS  ( Table  3.111 I 
In  experimenting  to  find  a  "good"  combination  of  I<(-  and  for  the  AAW  ship 
one  thing  is  clear.  <  ZA|.  since  the  AAW  ship  has  superior  AAW  capabili¬ 
ties.  This  implies  that  we  are  considering  a  subset  of  X.-  ,  such  that  I < 
ZA>  .  In  Experiment  3.1  of  the  single  ship,  single  sector  case,  we  found  that 

*  .2  and  *  1000  yielded  P^  *  .934  and  R  *  22.58;  reasonable  results. 

In  Table  3.11.1,  the  optimal  ship-  to-carrier  distance,  R,  increased  by  about 
5  miles.  This  is  expected  and  acceptable  since  the  ASW  threat  encourages  a 
ship  (AAW  or  ASW)  to  be  even  further  from  the  carrier  for  ASW  surveillance 
reasons.  This  is  a  result  of  the  *£  value  selected.  The  Pp  value  drops  because 
the  ship  is  faced  with  a  dual  threat. 

Based  on  the  results  of  Tables  3.111  we  choose  as  our  canonical  values  for  the 
AAW  ship: 

kA,*  3 

*  1 

IAi  *  1000 
and, 

!$,•  *  500 
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„v 


I 

I 


1 


’•  •] 
.  -< 


3.8  AN  ASW  SHIP  -  EXPERIMENTAL  RESULTS  (DUAL  CAPABILITIES) 


As  in  Section  3.7,  we  wish  to  maximize  the  probability  of  survival  of  the  car¬ 
rier,  under  an  AAW  and  an  ASW  threat.  Now,  however,  the  carrier  is  being 
defended  by  am  ASW,  ship.  We  impose  criteria  for  optimality:  the  optimal  proba¬ 
bility  of  survival,  must  exceed  .8  and  optimal  ship-to-carrier  distance,  R, 
must  be  approximately  50.  We  wish  to  characterize  an  ASW  ship  by  finding 
canonical  values  for  XA;  ,  X^j  ,  and  1 that  satisfy  these  criteria. 

3.8.1  SUITABLE  KAl  ,  K*;  ,  2^1  AND  Ski  VALUES 

We  consider  Eq’ns  (3.3)  and  (3.7)  which  represent  the  ASW  ship's  ability  to 
Xill  AAW  and  ASW  targets  respectively. 


Pv  =  Prv  • 


(“d\  > 


where , 


Pr  »  1  -  e 

A ; 


-(ka; c; > 


,  X  >  0 


d*.  *  xjsintf  -  yt- cos0 

1^1  is  the  AAW  range  factor  and, 


(3.3) 


Chapter  3  Single  Ship  -  Single  Sector  and  Dual  Capabilities  Case  Studies 


6  is  the  angle  of  approach  of  the  AAW  target. 

Also, 

(-d*1  /Z<.  ) 

Pft.  *  P e  '  '  (3.7) 

where , 

r  -  ) 

Pr  .  ■  1  -  e  1  ,  k  >  0 

d*  *  x-'Sin$  -  yt-cos $ 

£g>  is  the  ASW  range  factor  and, 
d  is  the  angle  of  approach  of  the  ASW  target. 

What  are  suitable  values  for  k.- ,  k_- ,  !*•  and  £-•  ? 

n*  ?i  M  •  Pi 

EXPERIMENT  3.12 

Purpose:  To  find  values  for  k^  ,  k5j  ,  Z^;  ,  and  Z^-  ,  such  that  the  final 
probability  of  survival  of  the  carrier,  P^.,  is  greater  than  .8 
and  R,  the  optimal  ship-to-carrier  distance  is  about  50. 

The  AAW  and  ASW  sectors  extend  from  0-45  and  both  have  a  single 
target  threat. 

Table  3.120  shows  the  values  of  the  fixed  parameters  in  Experiment  3.12. 
Table  3.121  shows  the  effect  of  varying  k^t-  . 
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Number 

Type 

Number 

22! 

Sectors 

Number 

tsi 

of 

of 

of 

of 

of 

Ships 

Ship 

Sectors 

Sectors 

Targets 

1 

ASW 

2 

1  AAW 

0-45 

1  AAW 

8000 

1  ASW 

0-45 

1  ASW 

Table  3.120  Fixed  Parameters  in  Exp't  3.12:  1  Ship,  Dual  Capabilities 


*a; 

800 

800 

850 

.05 

.06 

.06 

Surv.  Prob.  P 

.776 

.807 

.813 

Ship  Range  R 

47.75 

44.63 

45.17 

Ship  Angle  f 

22.  ? 

2? 

22.5 

Table  3.121  Effect  of  Varying  IAi  and  kA . 
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DISCUSSION  OF  RESULTS  (Table  3.121 ) 

We  know  that  as  k^.  increases,  increases.  Also  as  1^.  increases 

increases.  By  introducing  an  AAW  threat,  we  expect  the  optimal  ship-to-carrier 
distance,  R,  to  decrease.  Most  importantly,  we  are  restricted  to  a  subset  of 
ZA(-  values.  We  know  that  ZA;  must  be  less  than  8000  since  an  ASW  ship  has 
superior  ASW  capabilities.  Also,  ZA;  must  be  less  than  1000  =  Z*;  of  the  AAW 
ship.  If  not,  then  why  bother  having  AAW  ships  if  the  ASW  ships  can  be  more 
effective  against  AAW  targets? 

Base  on  the  results  of  Table  3.121  we  select  as  our  canonical  values  for  the 
ASW  ship: 

a  *2 

ks;  *  -1 

Z Ai  3  850 

and, 

Z5;  =  8000 
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We  now  have  established  characteristic  k^;  ,  k^  ,  S^-  ,  and  I values  for 
both  our  AAW  and  ASW  ships.  These  canonical  ships  will  be  used  in  Chapter  4  to 
investigate  the  behaviour  of  the  model  under  different  scenarios  when  two  or 
more  ships  are  protecting  the  carrier. 

3.9  CONCLUDING  REMARKS 

In  this  chapter,  we  performed  several  numerical  experiments.  The  purpose  of 
these  experiments  was  to  observe  the  effect  on  the  carrier  survival  probabili¬ 
ty  and  on  the  optimum  location  of  the  single  AAW  (or  ASW)  ship,  under  differ¬ 
ent  scenarios. 
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CHAPTER  4  MULTIPLE  SHIP  CASE  STUDIES  WITH  AAW  AND  ASW  THREAT  SECTORS 

4.0  SUMMARY 

We  venture  into  the  realm  of  two  ships.  We  investigate  how  the  two  ships  work 
in  harmony,  in  order  to  realize  the  maximum  survival  probability  of  the  carri¬ 
er.  We  demand  that  the  model  trek  into  hazardous  territory;  pure  threats,  mul¬ 
tiple  targets,  changes  in  the  sector,  tests  for  symmetry.  We  further  impose 
the  additional  burden  of  three  ships.  But  the  model  successfully  glides 
through  these  tests,  and  emerges  undaunted. 

4.1  PROBLEM  STATEMENT 

The  purpose  of  the  problem  is  to  maximize  the  probability  of  survival  of  a 
carrier  given  air  and/or  sub  threats.  The  threats  are  uniformly  distributed 
over  known  sectors.  An  AAW  and  an  ASW  ship  are  defending  the  carrier.  The  sol¬ 
ution  to  the  problem  is  the  optimum  locations  of  the  AAW  and  the  ASW  ships 
with  respect  to  the  carrier;  those  locations  that  maximize  the  probability  of 
survival  of  the  carrier. 
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The  purpose  of  this  chapter  is  to  study  various  scenarios  of  the  two  and  three 
ship  problems  in  order  to  see  how  the  ships  interact  to  best  protect  the  car¬ 
rier. 

4.2  MATHEMATICAL  FORMULATION ;  TWO  SHIPS 

The  carrier  is  positioned  at  the  origin.  The  AAW  ship  and  the  ASW  ship  have 
coordinates  (x|fyf)  and  respectively,  which  represent  their  locations 
relative  to  the  carrier.  In  polar  coordinates  (R, ,ft  )  denotes  the  optimal 
location  of  the  AAW  ship,  while  (Rz,fx)  denotes  the  optimal  location  of  the 
ASW  ship. 

4.2.1  OPTIMUM  INTERCEPTION  DISTANCE 

We  defined  the  (optimum)  interception  distance. to  be  the  shortest  distance 
between  ship  i  and  the  target  trajectory.  This  is  the  distance  at  which  the 
ship  attempts  to  kill  the  target.  In  Section  2.2.1,  we  said  that  d*- was  the 
intercept  distance  and  that 

d*.  =  x;sin0  -  yrcostf  (4.1) 
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where, 

(x’,yj)  are  the  coordinates  of  ship  i,  and 
$  is  the  angle  of  approach  of  the  target. 

4.2.2  THE  EFFECT  OF  THE  SHIPS'  DISTANCE  TO  THE  CARRIER 

In  Section  3.6.2  we  determined  that  the  ships  were  required  to  keep  a  certain 
distance  from  the  carrier.  We  introduced  the  functions  Pr  . ,  and  Pr  . ,  for 
every  ship  i,  to  capture  this  constraint.  From  Eq'ns  (3.4)  and  (3.6)  we  have 
that 


-«c  r  -  ) 

%,  *  1  -  «  ’  *  »  0  (4‘2> 


where, 

rs - ~ 

rj  *  Vx-(  +  y;  is  the  distance  of  ship  i  to  the  carrier,  and 

« 

is  a  parameter  associated  with  each  cruiser  and  destroyer. 


We  also  have 
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-(k*  ri  } 

P,  =  1  -  e  ,  lc  >  0 


(4.3) 


where , 


rj  *»x-(  +  y.  is  the  distance  of  ship  i  to  the  carrier, 
X, .  is  a  parameter. 


4.2.3  PROBABILITY  OF  INTERCEPTION 


As  we  explained  in  Section  2. 2. 2. 2,  each  ship  has  two  functions  which  fully 
describe  its  ability  to  intercept  targets.  P^.  of  Eq'n(4.4)  reflects  ship  i's 
ability  to  intercept  an  AAW  target. 


P^.  ,  the  PROBABILITY (SHIP  I  KILLS  AN  AAW  TARGET),  is 


("d»  /I  A-  ) 


PA-  3  prA  • 


where. 


(4.4) 


PrA;  -  1  ”  e 


-<*  ri  > 


'  v° 


d»  ■  (xj sin0  -  y^-cosfl), 

£..  is  the  AAW  range  factor,  and 

Ai 

6  is  the  angle  of  approach  of  the  AAW  target. 


(4.2) 


(4.1) 
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which  completely 


Similarly,  for  every  ship  i  we  have  a  function,  P<.  , 
describes  its  ability  to  intercept  ASW  targets. 


5; 

where , 


p«-  =  pr<-  e 


(-d*.  /I.) 

1  Ji 


-(**•  r;  } 

P  =  1-  e  '  ,  k  >  0 

S,  5. 


d*  *  x‘ sin0  -  y-cosd, 

i  1  • 

Z3;  is  the  ASW  range  factor,  and 
0  is  the  approach  of  the  ASW  target. 


(4.5) 


(4.3) 

(4.1) 


4.2.4  PROBABILITY  DISTRIBUTIONS  OF  THE  TARGETS 


Both  the  AAW  and  the  ASW  targets  are  uniformly  distributed  over  known  sectors. 
Let  (a,b)  represent  such  a  sector.  Let  f(9)  be  the  probability  density  func¬ 
tion  of  an  AAW  or  an  ASW  target.  Then  f(0)  is  given  by 

f(0)  •  l/(b  -  a)  a  <  9  *  b  (4.6) 

*  0  otherwise 
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4.2.5  SURVIVAL  PROBABILITY  OF  THE  CARRIER 


The  probability  that  the  carrier  survives  is  the  probability  that  at  least  one 
ship  kills  the  AAW  target  in  the  AAW  sector  (a,b)  and  at  least  one  ship  kills 
the  ASW  target  in  the  ASW  sector  (c,d). 


Prob (carrier  survives) 


(l  -  Pd-P*.  ))  id 


x 


(l 


na-p, 


i* 


s; 


))  d* 


(4.7) 


4.3  AM  AAW  AND  AN  ASW  SHIP  -  NUMERICAL  RESULTS 


We  wish  to  test  the  sensitivity  of  the  model  and  observe  how  the  two  ships 
work  in  conjunction  with  each  other  in  order  to  maximize  the  probability  of 
survival  of  the  carrier,  under  various  scenarios.  We  test  for  symmetry  and 
consider  the  result  of  a  pure  (AAW  or  ASW  but  not  both)  threat.  We  examine  the 
effects  of  multiple  targets  and  different  sectors  on  the  objective  function 
and  on  the  optimal  positions  of  tne  ships.  We  also  study  the  limiting  cases. 

We  use  the  following  notation  throughout  the  entire  chapter. 
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HOT  XT  I  OH 

The  maximum  probability  of  carrier  survival  =  . 

The  optimal  position  of  the  AAW  ship  from  the  carrier,  in  polar  coordinates, 
is:  (Rj  ,  )  • 

The  optimal  position  of  the  ASW  ship  from  the  carrier,  in  polar  coordinates, 
is :  (R^ ,  )  • 


4.3.1  SUITABLE  K M\  .  K<C  ,  Za\  .  AMP  2-s;  VALUES 

In  Experiments  3.11  and  3.12  of  Chapter  3,  we  discovered  appropriate  parameter 
values  to  represent  standard  AAW  and  ASW  ships  respectively.  They  are  shown 
in  Tables  4.0a  and  4.0b. 

These  canonical  values  will  be  used  throughout  this  chapter,  unless  otherwise 
specified. 

Note  that  the  ASW  ship  has  excellent  AAW  capability  (almost  as  good  as  the  AAW 
ship).  However,  the  AAW  ship  has  poor  ASW  capabilities  in  relation  to  those  of 
the  ASW  ship. 
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4.3.2  SINGLE  SHIP  -  SINGLE  THREAT  CASE 

We  return  for  a  moment  to  the  single  ship  -  single  threat  case  of  Chapter  3, 
in  order  that  we  may  use  it  as  a  point  of  comparison.  The  survival  probability 
of  the  carrier  in  Experiments  4.1  and  4.2  is  given  by  Eq'ns  (3.5)  and  (3.7) 
for  AAW  and  ASW  threats  respectively.  In  Experiment  4.1,  we  seek  the  values  of 
^  and  (R,  ,^,)  when  a  single  AAW  ship  is  required  to  defend  the  carrier 
against  an  AAW  threat  in  a  sector  extending  from  (o,  60).  Experiment  4.2  is 
analogous,  with  the  AAW  ship  replaced  by  an  ASW  ship. 

EXPERIMENT  4.1 

Purpose:  Given  the  canonical  k*,  and  ZA(  values  for  the  AAW  ship 
(Table  4.0a)  and  a  (0,60)  sector  with  a  single  AAW  target, 
to  determine  Pp  and  (R,  ,  ft  ) . 

Table  4.10  shows  the  value  of  the  fixed  parameters  in  Experiment  4.1. 

Table  4.11  shows  the  values  of  Pf  and  (R,,^  )• 

DISCUSSION  OF  RESULTS  (Table  4.W 

Because  of  the  uniform  disribution  of  the  target,  ft  =  30;  R(»  20.4  and 
Pp  *  .949  are  both  satisfactory. 
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Number 

Type 

Number 

Type  Sector 

Number 

of 

of 

of 

of 

of 

Ships 

Ship 

Sectors 

Sector 

Targets 

1 

AAW 

1 

AAW  0-60 

1 

Table  4.10  Fixed  Parameters  in  Fxp't  4.1:  1  Ship 


Surv.  Prob. 

pf 

.949 

Ship  Range 

R, 

20.40 

Ship  Angle 

t, 

30 

Table  4.11  Survival  Probability  and  AM  Ship  Location 
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EXPERIMENT  4.2 

Purpose:  Given  the  canonical  values  of  and  for  an  ASW  ship 

(Table  4.0b)  and  a  (0,60)  sector  with  a  single  ASW  threat, 
to  determine  the  values  of  Pp  and  (Rt*fz). 

Table  4.20  shows  the  values  of  the  fixed  parameters  in  Experiment  4.2. 

Table  4.21  shows  the  values  of  Pp  and 

DISCUSSION  OF  RESULTS  f Table  4.21) 

» 

Because  of  the  uniform  distribution  of  the  target,  fx-  30;  Rj=  46.31  and 
Pp  *  .968  are  both  satisfactory. 


4.3.3  SYMMETRY 


It  is  worthwhile  and  interesting  to  discover  what  role  symmetry  plays  in  the 
model.  The  next  experiment  uses  different  initial,  conditions  to  determine  how 
the  optimal  ship-to-carrier  distances  R,  and  R^,  and  their  associated  angles 
f,  and  +2*  are  affected. 
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Humber 

Type 

Number 

Type 

Sector  Number 

of 

of 

of 

of 

of 

Ships 

Ship 

Sectors 

Sector 

Targets 

1 

ASW 

1 

ASW 

0-60  1 

Table  4.20  Fixed  Parameters  in  Exp't  4.2:  1  Ship 


Surv.  Prob.  Pf 

.968 

Ship  Range  Rv 

46.41 

Ship  Angle  f 

30 

Table  A.21  Survival  Probability  and  ASH  Ship  Location 
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EXPERIMENT  4.3 

Purpose:  To  identify  under  what  conditions  we  may  expect 

symmetry.  We  use  the  parameter  values  in  Tables  4.0a  and  4.0b, 
and  identical  (0,60)  AAW  and  ASW  threat  sectors,  each  with  a 
single  target. 

Table  4.30  shows  the  values  of  the  fixed  parameters  in  Experiment  4.3. 

Table  4.31  shows  P^  ,  (R,,^()  and  (R*,?*.)  as  a  result  of  different 
initial  conditions.  (The  initial  conditions  are  in  Cartesian 
coordinates. ) 

DISCUSSION  OF  RESULTS  ( Table  4.31/ 

o 

As  we  would  expect,  there  is  symmetry  about  the  30  line.  See  Figures  4. la  and 
4.lb. 

This  is  apparent  when  we  compare  the  first  two  cases  of  Table  4.3i.  In  the 

o 

first  case,  the  AAW  ship  is  below  the  30  line  with  initial  coordinates  (20,5), 

o 

while  the  ASW  ship  is  above  the  30  line  with  initial  coordinates  (30,40).  The 

optimal  locations  reflect  a  small  change  in  positon,  and  =  15.02  and 

o 

fz=  45.3  indicates  that  the  AAW  ship  is  still  below  the  30  line  and  that  the 
ASW  ship  is  still  above  the  30°line.  R, =  24.51  and  Rx=  38.85  are  both  desira¬ 
ble  AAW  and  ASW  locations  respectively. 
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Number 

Type 

Number 

Type 

Sectors 

Number 

of 

of 

of 

of 

of 

Ships 

Ships 

Sectors 

Sectors 

Targets 

2 

1  AAW 

2 

1  AAW 

0-60 

1  AAW 

1  ASW 

1  ASW 

0-60 

1  ASW 

Table  4.30  Fixed  Parameters  in  Exp't  4.3:  2  Ships 


I nit.  Cond'n 
(Ship  Position) 
Surv.  Prob.  Pp 
(R,  ) 


Case  1 

(20,5)  AAW 
(30,40)  ASW 
.986 

(24.51,15.02) 

(38.85,45.30) 


Case  2 
(50,50)  AAW 

(10,5)  ASW 
.986 

(24.51,44.98) 

(38.85,14.70) 


Case  3 

(20,5)  AAW 
(40,10)  ASW 
.986 

(24.51,44.98) 

(38.85,14.70) 


Table  4.31  Survival  Probability  and  Ship  Locations 
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dentical  Threat  Sectors  (o  ,60) 
ion  of  AAW  Ship:  (24.51,  15. 02°. 
ion  of  ASW  Ship:  (38.85,  45.3°) 


One  AAW,  One  ASW 


of  Symmetry  (along  with  Figure  4.1b) 


.ASVs /  SHIP 


0  ZD 


40 


'AUtfS 


~.tical  Threat 
i  of  AAW  Ship: 
i  of  ASW  Ship: 


Sectors 
'24. 51, 
(38.85, 


(0  ,60): 

o 

44.98  ) 
o 

14.7  ) 


One  AAW,  One  ASW 


Example  of  Symmetry  (along  with  Figure  4.1  a  I 
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In  the  second  case,  we  use  an  initial  condition  of  (50,50)  for  the  AAW  ship, 

and  (10,5)  for  the  ASW  ship.  Now  the  AAW  ship  is  on  the  45°line,  and  the  ASW 

o  o 

ship  is  below  the  30  line.  The  final  result  of  (24.51,44.98)  for  the  AAW  ship, 

o 

and  (38.85,14.7)  for  the  ASW  ship  demonstrates  that  the  image  was  produced 

e 

across  the  30  line,  and  therefore  we  have  symmetry.  It  is  also  interesting  to 
note  that,  given  the  initial  conditions,  the  model  brought  in  the  AAW  ship  and 
pulled  out  the  ASW  ship  as  was  desired.  The  final  objective  function  value  ^ 
is,  of  course,  the  same  in  both  cases. 

4.3.4  PURE  THREATS  -  TWO  SHIPS 

We  now  turn  to  the  situation  in  which  only  one  type  of  threat  exists  -  AAW  or 
ASW,  but  not  both.  The  carrier  is  still,  however,  being  defended  by  both  the 

AAW  and  the  ASW  ship.  This  scenario  is  of  interest  since  it  may  shed  light  on 

the  behaviour  of  the  model  under  multiple  AAW  (or  ASW,  but  not  both)  threats. 

For  example,  if  the  number  of  AAW  threats  is  increased  to  100,  and  there  is  a 

single  ASW  threat,  do  the  optimal  ship-to-carrier  locations  (R, )  and 
converge  to  those  of  the  situation  in  which  there  is  only  a  pure  AAW 
threat?  In  order  to  answer  this  question,  we  need  the  results  of  (R,,^,  )  and 

(Rx» ti)  for  tfie  pure  threat;  hence,  the  following  two  experiments.  It  is  suf¬ 
ficient  to  run  the  experiments  with  only  one  target,  since  we  saw  in  Claim  3.1 
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that  the  optimal  ship-to-carrier  distance  is  independent  of  the  number  of  tar¬ 
gets. 

EXPERIMENT  4.4 

Purpose:  To  obtain  the  values  of  (R,,^t)  and  (R z'fi)  for  a  Pure  AAW 
threat.  There  is  only  one  (0,60)  sector  with  a  single 
AAW  threat. 

Table  4.40  shows  the  values  of  the  fixed  parameters  of  Experiment  4.4. 

Table  4.41  contains  the  results  of  Experiment  4.4;  Pf  ,  (R(,?()  and 

^2/^2. )  • 

DISCUSSION  OF  RESULTS  < Table  4.41) 

Because  both  the  AAW  and  the  ASW  ship  have  excellent  AAW  capabilities 
(XA|  ■  1000,1^=  850),  the  survival  probability,  P^»  ,  is  very  high. 

R,  =  22.10  and  21.49  both  reflect  typical  AAW  ship  locations,  i.e. 

o  o 

locations  which  are  close  to  the  carrier,  =,  15.21  and  f 44.91  reveal 

o 

that  the  ships  are  well  positioned  about  the  30  line.  We  know  that,  given  dif¬ 
ferent  initial  conditions,  we  could  obtain  symmetric  optimal 
locations, (results  of  Experiment  4.3).  See  Figure  4.2. 
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Number 

Type 

Number 

Type  Sector 

Number 

of 

of 

of 

of 

of 

Ships 

Ships 

Sectors 

Sector 

Targets 

2 

1  AAW 

1 

AAW  0-60 

l 

1  ASW 

Table  4.40  Fixed  Parameters  in  Exp't  4.4 :  2  Ships 


Surv.  Prob. 


(Ri ) 

(1 H't-L ) 


.997 

(22.10,15.21) 

(21.49,44.91) 


Table  4.41  Survival  Probability  and  Ship  Locations :  Pure  AM  Threat 
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EXPERIMENT  4.5 

Purpose:  To  obtain  the  values  of  (R ( t )  and  (1^#^)  for  a  pure  ASW 
threat  in  a  (0,60)  sector. 


Table  4.50  shows  the  values  of  the  fixed  parameters  in  Experiment  4.5. 
Table  4.51  contains  the  results  of  Experiment  4.5;  P^-  ,  (R,,^),  and 
* 


DISCUSSION  OF  RESULTS  (Table  4.511 

The  ASW  ship  has  an  extremely  high  ASW  capability  and  this  is  reflected  in  a 
high  survival  probability.  Its  optimal  ship-to-carrier  distance,  R^,  is  typi¬ 
cal  of  that  of  an  ASW  ship  in  that  it  is  relatively  far  from  the  carrier.  In 
contrast,  the  AAW  ship  has  a  fairly  poor  ASW  capability,  (I5  *  500)  and  its 
optimal  ship-to-carrier  distance  represents  that  of  an  AAW  ship,  suggesting 

that  the  large  survival  probability  is  largely  due  to  the  ASW  ship.  Both  snips 

o 

are,  nevertheless,  well  positioned  about  the  30  line.  See  Figure  4.3. 
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Number 

Type 

Number 

Type 

Sector  Number 

of 

of 

of 

of 

of 

Ships 

Ships 

Sectors 

Sector 

Targets 

2 

1  AAW 

1 

1  ASW 

0-60  1 

1  ASW 


Table  4. SO  Fixed  Paraataters  in  Exp't  4.5:  2  Ships 


Surv.  Prob.  P„ 


(R.  *tt  > 


.934 

(29.40,45.98) 

(49.10,16.92) 


Table  4.51  Survival  Probability  and  Ship  Locations:  Pure  ASH  Threat 
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Legend 

1) 

Location 

of 

AAW  Ship: 

(29.4 

2) 

Location 

of 

ASW  Ship: 

(49.1 

3) 

Threat  Sector  (o  ,  60 

) 

45.98°) 


16.92°) 


ti 


Figure  4.3 


Locations  of  Ships  for  Pure  ASH  Threat 


y. 
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4.3.5  MULTIPLE  TARGETS,  LIMITING  CASES,  AND  A  CHANGE  IN  SECTOR 

In  this  section,  we  examine  the  consequences  of  multiple  targets,  probe  the 
limiting  cases  and  study  the  results  of  a  change  in  sector.  We  want  to 
observe  how  the  two  ships  interact  to  produce  the  maximum  probability  of  sur¬ 
vival  of  the  carrier  under  these  types  of  scenarios. 

EXPERIMENT  4.6 

Purpose:  To  obtain  the  optimal  locations,  i.e.  the  values  of  (R,  ,f.  )  and 
( Rj, f  j)  for  different  numbers  of  AAW  threats,  while  the  number 
of  ASW  threats  is  held  constant  at  one.  The  AAW  threat  sector 
and  the  ASW  threat  sector  both  extend  from  (0,60). 

Table  4.60  shows  the  values  of  the  fixed  parameters  of  Experiment  4.6. 

Table  4.61  shows  the  effect  on  Pp  ,  (R,  ,f%  )  and  (Rz,^2)  of  changing  the 
AAW  threat  concentration. 

DISCUSSION  OF  RESULTS  (Table  4.61 ) 

The  results  of  the  experiment  show  that  as  the  number  of  AAW  targets  increase, 
R,  and  Rjl  increase,  as  do  fy  and  This  information  is  particularly  inter¬ 
esting  when  we  compare  it  to  the  results  of  Experiment  4.4.  In  the  limiting 
case,  when  the  number  of  AAW  targets  is  "large",  the  single  ASW  target  is  of 
little  importance,  as  the  two  ships  virtually  confine  their  intentions  to  the 
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Number 

of 

Ships 

Type 

Of 

Ships 

Number 

of 

Sectors 

Type 

of 

Sectors 

Sectors  Number 
of 

ASW  Targets 

2 

1  AAW 

2 

1  AAW 

0-60 

l 

1  ASW 

1  ASW 

0-60 

Table  4.60  Fixed  Parameters  in  Exp't  4.6: 

2  Ships 

Number  of  AAW  Targets 

Pf 

<v* > 

1 

.985 

(24.51,44.98) 

(38.80,14.70) 

2 

.980 

(24.00,45.00) 

(35.00,14.40) 

3 

.975 

(23.80,45.10) 

(32.50,14.38) 

6 

.962 

(23.40,45.10) 

(28.60,14.41) 

12 

.940 

(23.00,45.00) 

(25.50,14.50) 

25 

.899 

(22.70,44.97) 

(23.60,14.73) 

50 

.827 

(22.48,44.90) 

(22.60,14.87) 

Table  4. SI  Effect  of  Varying  the  Number  of  AM  Targets 
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massive  AAW  threat.  The  optimal  survival  probability,  ,  is  almost  com¬ 
pletely  determined  by  the  AAW  capability  of  each  ship,  and  .  We  see 

by  the  AAW  capability  of  each  ship,  and  .  we  see  that  (Rj.^and 

(Rj,^)  of  Experiment  4.6  converge  to  (Rj  ,f%  )  and  (Rz»?x)  of  Experiment  4.4, 
respectively,  in  the  limit. 

Note  that  in  this  experiment,  the  major  effect  of  increasing  the  number  of  AAW 
targets  is  in  the  location  of  the  ASW  ship.  For  increasing  number  of  targets, 
the  ASW  ship  is  being  pulled  more  and  more  toward  the  carrier,  very  minor 
changes  occur  in  the  location  of  the  AAW  ship.  This  is  due  to  the  excellent 
AAW  capability  of  the  ASW  ship. 

EXPERIMENT  4.7 

Purposes  To  obtain  the  values  of  (R,,^,)  and  (R uti)  under  different 
ASW  target  densities,  when  the  number  of  AAW  targets  is  held 
constant  at  one.  The  AAW  and  the  ASW  threat  sectors  both  extend 
from  (0,60). 

Table  4.70  contains  the  values  of  the  fixed  parameters  of  Experiment  4.7. 

Table  4.71  shows  the  effect  on  ,  (R( ,fx  )  and  (R^,  H)  of  changing  the 
number  of  ASW  targets. 
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Number 

of 

Ships 

Type 

of 

Ships 

Number 

of 

Sectors 

Type 

of 

Sectors 

Sectors 

Number 

of 

AAW  Targets 

2 

1  AAW 

2 

1  AAW 

0-60 

1 

1  ASW 

1  ASW 

0-60 

Table  4.70  Fixed  Parameters  in  Exp't  4.7:  2  Ships 


Number  of  ASW  Targets 


1 

.985 

2 

.979 

3 

.972 

6 

.953 

12 

.917 

25 

.844 

50 

.721 

100 

.527 

(R| »f,  ) 

(V'i.) 

(24.51,44.98) 

(38.85,14.70) 

(25.26,44.97) 

(42.11,15.20) 

(25.84,45.02) 

(43.62,15.51) 

(26.85,45.18) 

(45.65,16.04) 

(27.70,45.40) 

(47.00,16.40) 

(28.40,45.60) 

(47.90,16.60) 

(28.90,45.70) 

(48.50,16.70) 

(29.10,45.80) 

(48.70,16.80) 

Table  4.71  Effect  of  Varying  the  Number  of  ASM  Targets 
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DISCUSSION  OF  RESULTS  (Table  4.71) 

It  is  no  surprise  that  we  will  draw  the  same  conclusions  here  chat  we  did  for 
Experiment  4.6.  Again,  (R, ,?,)  and  of  this  experiment  converge  to 

(R, ,^t)  and  of  Experiment  4.5,  respectively,  in  the  limit,  as  the  num¬ 

ber  of  ASW  targets  gets  large.  This  is  so  because  as  the  ASW  threat  concen¬ 
tration  increases,  the  AAW  threat  becomes  negligible,  and  the  two  ships  are 
required  to  flex  their  ASW  capability  muscles. 

In  this  experiment,  as  the  number  of  ASW  targets  increases,  both  ships  are 
moved  away  from  the  carrier  in  a  coordinated  manner.  Even  though  the  AAW  ship 
has  limited  ASW  capabilities,  they  are  still  being  effectively  utilized. 

Next,  we  see  what  happens  when. the  ASW  sector  is  shifted. 

EXPERIMENT  4.8 

Purpose:  To  see  the  effect  on  Pf  ,  (R,,f,)  and  of  shifting  the 

ASW  sector.  The  AAW  sector  does  not  change  and  extends  from 
(0,60).  Both  the  AAW  sector  and  the  ASW  sector  have  a  single 
target  threat. 

Table  4.80  shows  the  values  of  the  fixed  parameters  of  Experiment  4.8. 

Table  4.81  shows  the  effect  on  Pf  ,  (R,,f,  )  and  (Rj_, f j_)  of  changing  the 
ASW  sector. 
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Number 

of 

Ships 

of 

Ships 

Number 

of 

Sectors 

Type 

of 

Sectors 

AAW  Sector 

Number 

of 

Targets 

2 

1  AAW 

2 

1  AAW 

0-60 

1  AAW 

1  ASW 

1  ASW 

1  ASW 

Table  4.80  Fixed  Parameters  in  Exp't  4.8:  2  Ships 


ASW  Sector 

Pf 

(R|,f.  ) 

(RV?l) 

0-60 

.985 

(24.51,15.02) 

(38.80,45.30) 

10-70 

.984 

(23.41,16.84) 

(40.81,47.69) 

20-80 

.981 

(21.91,18.26) 

(40.61,50.68) 

30-90 

.975 

(20.60,19.56) 

(39.29,54.51) 

Table  4.81  Effect  of  Changing  the  ASH  Threat  Sector 
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DISCUSSION  OF  RESULTS  (Table  4.81) 

As  the  ASW  sector  was  shifted  from  (0,60)  to  (30,90),  the  AAW  ship  is  being 
pulled  in  toward  the  carrier  and  moves  more  and  more  toward  the  center  of  the 
AAW  sector,  to  provide  better  AAW  coverage.  The  range  of  the  ASW  ship  does  not 
change  much;  its  angle  changes  significantly  moving  it  toward  the  center  of 
the  ASW  threat  sector.  However,  the  ASW  ship  still  contributes  to  the  AAW 
defense.  A  decrease  in  the  optimal  survival  probability,  P^.  ,  as  the  threat 
sector  area  expands,  also  supports  the  theory  that  the  two  ships  cannot  defend 
the  carrier  as  well. 

What  happens  when  we  increase  the  number  of  AAW  targets  under  this  type  of 
scenario?  This  is  the  subject  of  the  next  experiment. 

EXPERIMENT  4.9 

Purpose;  To  determine  the  effect  on  (Rj  )  and  (R^,^,)  of  increasing 

the  number  of  AAW  targets,  while  keeping  the  ASW  target  constant 
at  one.  The  AAW  sector  extends  from  (0,60l  and  the  ASW  sector 
is  from  (30,90). 

Table  4.90  contains  the  values  of  the  fixed  parameters  of  Experiment  4.9. 

Table  4.91  demonstrates  the  effect  of  increasing  the  number  of  AAW 
targets. 
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i 

i 


Number 

Type 

Number 

Type 

Sectors 

Number 

of 

of 

of 

of 

of 

Ships 

Ships 

Sectors 

Sectors 

ASW  Targets 

2 

1  AAW 

2 

1  AAW 

0-60  AAW 

1 

1  ASW 

1  ASW 

30-90  ASW 

Table  4.90  Fixed  Parameters  in  Exp't  4.9:  2  Ships 


i 


Number  of  AAW  Targets  Pp 


1 

.975 

2 

.967 

4 

.955 

8 

.937 

16 

.907 

32 

.856 

(RrfS  ) 

vti 

(20.60,19.56) 

(39.29,54.51) 

(21.04,17.32) 

(36.38,50.55) 

(21.38,16.05) 

(33.12,48.01) 

(21.61,15.41) 

(29.65,46.53) 

(21.75,15.18) 

(26.56,45.71) 

(21.85,15.14) 

(24.32,45.28) 

Table  4.91  Effect  of  Varying  the  Number  of  AAI4  Targets 
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DISCUSSION  OF  RESULTS  (Table  4.91) 

As  the  AAW  threat  becomes  large,  the  ASW  threat  becomes  negligible.  The  opti¬ 
mal  positions  of  the  ships  are  dictated  essentially  by  the  existence  of  the 
AAW  threat,  and  we  notice  that  (R,,^,  )  and  of  this  experiment  converge 
to  (R, , ft  )  and  (R^,?^)  of  Experiment  4.5,  (pure  AAW  threat)  respectively,  in 
the  limit  as  the  number  of  AAW  targets  gets  large. 

Note  that  the  range  of  the  AAW  ship  does  not  change  significantly;  however, 
its  location  shifts  to  protect  the  lower  half  of  the  AAW  threat  sector.  The 
ASW  ship  is  pulled  close  to  the  carrier  and  moves  to  protect  the  upper  half  of 
the  AAW  threat  sector. 


4.3.6  THREE  SHIPS 


We  talce  a  glimpse  at  how  the  model  operates  when  confronted  with  three  ships 
protecting  the  carrier.  We  choose  to  observe  its  behaviour  with  two  AAW  ships 
and  one  ASW  ship.  In  the  succeeding  experiments,  we  consider  a  single  AAW  sec¬ 
tor  extending  from  (0,60)  and  a  single  ASW  sector  extending  from  (30,90).  We 
want  to  establish  how  the  three  ships  interact  in  order  to  produce  the  maximum 
survival  probability  of  the  carrier,  given  different  target  densities.  The  ASW 
ship  and  one  of  the  AAW  ships  have  the  parameter  values  of  Tables  4.0a  and 
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4.0b  respectively.  For  the  second  AAW  ship,  we  see  the  effect  of  varying  two 
of  its  parameters. 

The  pertinent  mathematical  equations  follow  directly  from  Section  2. 2. 3.1  Case 
4. 


NOTATION 

The  optimal  ship-to-carrier  location  of  the  second  AAW  ship  = 

EXPERIMENT  4.10 

Purpose:  To  see  the  effect  on  ,  (R,,ft),  (Rz,?z)  and  (R3,?5)  of 
increasing  the  number  of  ASW  targets,  when  the  second  AAW 
ship  has  essentially  no  ASW  capability. 

Table  4.100  shows  the  values  of  the  fixed  parameters  of  Experiment  4.10. 
Table  4.101  shows  the  result  of  increasing  the  ASW  target, 
concentration. 
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Number 

Type 

Number 

Type 

Sectors 

Number 

of 

of 

of 

of 

of 

Ships 

Ships 

Sectors 

Sectors 

Targets 

3 

2  AAW 

2 

1  AAW 

0-60  AAW 

1 

1  ASW 

1  ASW 

30-90  ASW 

-A*  -S5  §S5 

.2  0  1000  10 


Table  4.100  Fixed  Parameters  in  Exp't  4.10:  3  Ships 


ASW  Targets 

Pf 

'h  > 

l 

.992 

(25.29,43.33) 

(46.39,69.82) 

(25.13,11.42) 

10 

.937 

(28.79,43.97) 

(48.76,72.73) 

(24.98,11.67) 

20 

.880 

(29.09,43.99) 

(48.93,72.91) 

(24.98,11.68) 

Table  4.101  Effect  of  Varying  the  Number  of  ASH  Targets 
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DISCUSSION  OF  RESULTS  (Table  4.101) 

As  the  number  of  ASW  targets  increases  from  one  to  ten,  the  first  AAW  ship  is 
called  to  the  aid  of  the  ASW  ship;  it  is  moved  out  and  its  angle  is  increased, 
as  its  ASW  capability  becomes  more  useful.  The  second  AAW  ship,  due  to  its 
extremely  poor  ASW  performance,  is  hardly  moved  at  all.  As  the  number  of  ASW 
targets  is  yet  further  augmented  to  twenty,  we  witness  even  less  of  a  change 
in  optimal  positions,  suggesting  that  we  are  converging  to  the  best  locations. 
See  Figure  4.4. 

Next,  we  improve  the  ASW  capability  of  the  second  AAW  ship  . 

EXPERIMENT  4.11 

Purpose:  To  see  the  effect  on  ,  (R(,f,),  (R^'fi)  and 

of  increasing  the  number  of  ASW  targets,  when  the  second  AAW 
ship  has  a  good  ASW  capability,  (better  than  that  of  the  first 
AAW  ship) . 

Table  4.110  shows  the  values  of  the  fixed  parameters  of  Experiment  4.11. 

Table  4.111  shows  the  effect  on  P^  ,  (R, , ^  ),  (R^, tjJ  and  (Rj,^)  of 
increasing  the  ASW  target  density. 
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'X/fcWSHlP 

•ASW5UP 


2y9  AM  6*,p 

*  #2ypAAvvJ  5HlP 


zs 

50 

(25.29, 

43.33°) 

(1  Target) 

(46.39, 

o 

69.82  ) 

(1  Target) 

(25.13, 

11.42°) 

(1  Target) 

(29.09, 

43.99°) 

(20  Targets) 

(48.93, 

72.91° ) 

(20  Targets) 

* 


MIL£S 


M 
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Number 

of 

Ships 

Type 

of 

Ships 

Number 

of 

Sectors 

of 

Sectors 

Sectors 

Number 

of 

Targets 

"At,  «5s; t. 

3 

2  AAW 

2 

1  AAW 

0-60  AAW 

1  AAW 

.2  .1  1000  3000 

1  ASW 

1  ASW 

30-90  ASW 

Table  4.110  Fixed  Parameters  in  Exp't  4.11 :  3  Ships 


ASW  Targets 

Pf 

fl) 

1 

.997 

(21.08,13.10) 

(43.66,63.95) 

10 

.987 

(20.90,18.49) 

(49.14,74.39) 

20 

.980 

(20.59,22.50) 

(49.93,75.69) 

<«j  'V 

(25.29,35.47) 

(33.84,45.30) 

(35.66,46.87) 


Table  4.111  Effect  of  Varying  the  Number  of  ASH  Targets 
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DISCUSSION  OF  RESULTS  IT  able  4.111) 

Immediately  we  experience  the  effect  of  giving  the  second  AAW  ship  a  powerful 
ASW  capability.  Even  when  there  is  only  one  AAW  and  one  ASW  threat,  the  second 
AAW  ship  is  thrown  into  the  ASW  threat  sector  =  35.47)  emphasizing  its 
superior  ASW  capability.  The  first  AAW  ship  with  its  significantly  weaker  ASW 
capability,  is  brought  deeper  into  the  AAW  sector  (^(  =  13.10). 

It  is  useful  to  compare  the  results  of  Experiments  4.io  and  4.11  in  the  cases 
where  there  is  only  one  AAW  and  one  ASW  target.  In  Experiment  4.10,  the  first 
AAW  ship  acts  to  help  the  ASW  ship  because  the  second  AAW  ship  has  virtually 
no  ASW  capability.  In  Experiment  4.11,  the  roles  of  the  two  AAW  ships  are 
reversed,  because  the  second  AAW  ship  has  a  better  ASW  capability.  Also,  when 
the  number  of  ASW  targets  is  increased  from  one  to  twenty,  increases  sub¬ 
stantially,  demonstrating  that  the  second  AAW  ship  is  effectively  serving  as 
an  ASW  ship.  Its  angle,  f is  always  increasing,  bringing  it  more  towards  the 
ASW  sector. 

It  is  worthwhile  to  note  that,  in  Experiment  4.11,  the  optimal  survival  proba¬ 
bility,  P^.  ,  barely  suffers  as  the  ASW  targets  become  thicker.  It  decreases 
from  .997  (one  target)  to  .980  (twenty  targets).  On  the  other  hand,  the  lack 
of  ASW  capability  in  Experiment  4.10,  is  more  adequately  represented,  as  P^. 
plummets  from  .992  (one  target)  to  .880  (twenty  targets) . 
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EXPERIMENT  4.12 

Purpose:  To  see  the  effect  on  P f  ,  (R,  ),  (R^,^,)  and  (Rj, of 

augmenting  the  number  of  AAW  targets. 

Table  4.120  shows  the  values  of  the  fixed  parameters  of  Experiment  4.12. 

Table  4.121  shows  the  result  of  increasing  the  number  of  AAW  targets. 

DISCUSSION  OF  RESULTS  (Table  4.121) 

The  ASW  ship  has  an  excellent  AAW  capability.  As  soon  as  the  AAW  threat 
becomes  more  prominent,  the  ASW  ship  is  brought  down  towards  the  AAW  sector 
is  decreasing).  It  is  placed  closer  to  the  carrier,  talcing  on  the  optimal 
ship-to-carrier  distance,  Rz  ,  of  an  AAW  ship.  Because  all  three  ships  are 
equipped  with  good  AAW  capabilities,  the  optimal  survival  probability  remains 
excellent,  (P^  =  .990)  even  when  the  number  of  AAW  targets  is  increased  to 
twenty . 

How  good  are  the  results  when  there  are  three. threat  sectors?  As  a  final 
experiment,  we  consider  this  possibility. 
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Number 

of 

Type 

of 

Number 

of 

Type 

of 

Sectors 

Number 

of 

Ships 

Ships 

Sectors 

Sectors 

Targets 

3 

2  AAW 

2 

1  AAW 

0-60  AAW 

1  ASW 

.2  .1  1000 

3000 

1  ASW 

1  ASW 

30-90  ASW 

Table  4.120  Fixed  Parameters  in  Exp't  4.12:  3  Ships 


AAW  Targets 

Pf 

(Rj  tfx  ) 

(Rj.'  fl) 

(r3'^3) 

1 

.997 

(21.08,13.10) 

(43.66,63.95) 

(25.29.35.47) 

10 

.993 

(21.61,10.84) 

(35.59,53.05) 

(23.20,26.15) 

20 

.990 

(21.40,11.43) 

(32.09,51.98) 

(22.61,24.33) 

Table  4.121  Effect  of  Varying  the  Number  of  AW  Targets 
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EXPERIMENT  4.13 


Purpose:  To  see  the  effect  on  and  ( )  when 

the  number  of  ASW  targets  is  augmented.  The  ASW  sector  extends 
from  (60,90),  and  the  two  AAW  sectors  extend  from  (0,20)  and 
(30,50) . 

Table  4.130  shows  the  values  of  the  fixed  parameters  of  Experiment  4.13. 

Table  4.131  shows  the  results  of  Experiment  4.13. 

DISCUSSION  OF  RESULTS  f Table  4.131) 

For  a  single  AAW  and  ASW  threat,  the  model  successfully  places  an  AAW  ship  in 


each 

AAW 

sector, 

and 

the 

ASW  ship  in  the  ASW  sector 

.  It  also  wisely  places 

each 

ship 

roughly 

in 

the 

middle  of  its  sector.  The 

survival 

probability  is 

high; 

.996.  i 

When 

the 

ASW  target  threat  is  increased  to 

twenty,  all  the 

ships 

are 

shifted 

upwards 

toward  the  ASW  sector.  The 

optimal 

survival  proba- 

bility  decreases  (^.  =  .978)  . 


4.4  CONCLUDING  REMARKS 


In  this  chapter,  we  performed  numerical  experiments  to  investigate  how  two  and 
three  ships  worit  in  conjunction  with  one  another  in  order  to  maximize  the  sur¬ 
vival  probability  of  the  carrier. 
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Number 

Type 

Number 

Type 

Sectors 

Number 

of 

of 

of 

of 

of 

Ships 

Ships 

Sectors 

Sectors 

Targets 

3 

2  AAW 

3 

2  AAW 

0-20 

AAW 

1  AAW 

1  ASW 

1  ASW 

30-50 

AAW 

60-90 

ASW 

W  2*J  ^>4 

.2  .1  1000  3000 


Table  4.130  Fixed  parameters  in  Exp't  4.13:  3  Ships 


ASW  Targets 

Pf 

(R5'^S) 

1 

.996 

(24.65,  9.51) 

(48.30,74.06) 

(24.70,38.55) 

20 

.978 

(20.66,12.68) 

(57.27,77.78) 

(28.12,47.23) 

Table  4.131  Effect  of  Varying  the  Number  of  ASM  Targets 
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CHAPTER  5  CONCLUSIONS  AND  SUGGESTIONS  FOR  FUTURE  RESEARCH 


5.1  CONCLUSIONS 

A  simple  model  that  can  be  used  to  carry  out  quantitative  investigations  of 
the  optimal  location  of  multi-warfare  capable  platforms  to  carry  out  the  AAW 
and  ASW  functions  in  defense  of  a  carrier,  has  been  developed. 

Even  with  the  obvious  limitations  of  our  model,  the  results  of  the  research 
suggest  that  the  methodology  developed  can  be  used  in  conjunction  with  more 
realistic  models  of  AAW  and  ASW  capabilities;  in  this  case  it  is  possible  to 
develop  a  valuable  decision-aid. 

In  spite  of  its  simplicity,  the  model  captures  correctly  what  should  happen  in 
scenarios  in  which  the  coordination  of  multi-warfare  capable  platforms  is 
essential  for  superior  defense  against  simultaneous  AAW  and  ASW  threats.  This 
is  primarily  reflected  in  the  fact  that  the  optimal  ship  locations  with 
respect  to  the  carrier  change  in  a  reasonable  and  sensible  manner  as  a  func¬ 
tion  of: 

a)  the  location  and  number  of  AAW  and  ASW  threat  sectors 

b)  the  density  of  the  AAW  and  the  ASW  targets 

c)  the  individual  AAW  and  ASW  capabilities  of  each  of  the  ships 


Chapter  5  Conclusions  and  Suggestions  for  Future  Research 


PAGE  149 


The  nature  of  this  information  is  valuable  for  both  planning  and  tactical  sit¬ 
uations  involving  modern  Battle  Groups. 

The  model  is  a  complete  self-contained  unit;  it  is  mathematically  rigorous  and 
gives  consistent  results.  It  is  robust  in  the  sense  that  "small"  changes  in 
the  parameters  lead  to  "sensible"  changes  in  the  objective  function  and  the 
solution.  It  is  very  general  and  does  not  confine  us  to  a  narrow  stream  of 
information,  but  rather  opens  the  door  to  investigate  a  broad  cross-section  of 
scenarios.  We  are  free  to  select  any  combination  and  strength  of  AAW  and  ASW 
ships,  the  type  and  intensity  of  the  targets,  and  the  type,  location  and  num¬ 
ber  of  threat  sectors. 

A  special  advantage  of  the  model  is  that  it  is  user  -  friendly.  It  is 
extremely  easy  to  start  new  experiments,  and  the  results  are  clear.  Of  course, 
because  we  depend  on  numerical  approximations  via  the  optimization  routine 
(SUMSNO)  (5)  and  the  integration  routine  (D01BAF)  (6) ,  the  numerical  complexi¬ 
ty  of  our  model  is  predetermined. 


5.2  SUGGESTIONS  FOR  FUTURE  RESEARCH 
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The  intent  of  this  section  is  to  present  natural  modifications  and  exciting 
extensions  of  the  model.  We  propose  generalizing  the  sectors  to  include  not 
only  the  permissible  sectors,  but  also  combinations  of  non-permissible  and 
overlapping  multiple  sectors.  We  consider  the  possibility  of  diverse  probabil¬ 
ity  distributions  of  the  targets.  We  recommend  introducing  the  carrier  as  a 
ship  with  its  own  defense  capabilities.  Finally,  we  shift  gears  and  regard  the 
model  in  a  whole  other  dimension.  We  suggest  adapting  it  to  the  setting  of  a 
game,  and  create  the  conditions  under  which  it  can  be  used. 


5.2.1  PERMISSIBLE  AND  NON-PERMISSIBLE  AND  MULTIPLE  SECTORS 


In  Chapter  2,  we  developed  a  model  to  represent  the  survival  probability  of  a 
carrier,  given  air  and/or  sub  threats,  and  simple  models  of  the  probabilities. 
The  problem  was  formulated  in  two  dimensions,  with  the  carrier  at  the  origin. 
Our  objective  was  to  find  an  optimal  location  for  all  of  the  ships  defending 
the  carrier,  with  respect  to  the  carrier.  In  Section  2. 2. 3.1,  we  considered 
two  different  types  of  threat  sectors,  and  defined  the  notion  of  a  PERMISSIBLE 
sector  (DEF'N  2.9).  Permissible  sectors  are  ones  in  which  any  ship  wiil 
attempt  to  Kill  a  target  BEFORE  it  passes  through  the  origin. 

A  natural  way  to  extend  the  ideas  of  this  paper  would  be  to  incorporate  the 
possibility  of  non-permissible  sectors.  As  we  saw  in  Section  2. 2. 3.1,  a  nec- 
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essary  and  sufficient  condition  for  a  permissible  sector  (a,b)  is  that,  for  a 


given  ship  i  with  coordinates  (x^  ,y()  we  have 

d*.  =  x.'sin#  -  y-cos#  for  all  $  ,  a  ^  9  s  b 

i  1  i 

We  now  wish  to  relax  this  condition,  and  use  the  more  general  definition  of 
d*  ,  found  in  DEFINITION  2.8. 

i 

If  we  consider  the  most  general  scenario,  like  that  of  Case  4  at  the  end  of 
Chapter  2,  we  need  to  know  how  the  survival  probability  equation  is  affected 
by  the  existence  of  non-permissible  sectors.  First,  at  each  iteration  we  must 
determine  the  value  of  d«  .  In  DEFINITION  2.8,  we  stated  that 

d*.  =  x;sintf  -  y;  cosfl  ft-S  ft  $ 

=  rj  otherwise 


where , 

ft  is  defined  in  DEFINITION  2.6,  and 

and  are  defined  in  DEFINITION  2.7. 

It  is  easy  to  incorporate  this  information  into  the  computer  program,  it  suf¬ 
fices  to  declare  and  define  ft  ,  f, ,  and  and  to  add,  wherever  necessary,  a 
control  statement  of  the  form: 
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THEN  d*  =  x; sin#  -  y- cos# 

ELSE  d«.  =  r: 

i  i 

The  use  of  non-perraissible  sectors,  and  hence  the  inclusion  of  the  generalized 
d*  ,  leads  us  to  yet  a  further  extension  of  the  model.  We  are  now  free  to 
admit  multiple  overlapping  sectors.  The  ultimate  generalization  of  the  model 
in  this  way  is: 

Case  5:  H  Ships,  N{  Independent  and  Identically  Distributed  AAId  Targets 
in  AMd  Sector  i,  i  -  1,6,  Qj  Independent  and  Identically 
Distributed  ASId  Targets  in  ASId  Sector  j,  j  -  1,H. 

Then  the  survival  probability  becomes: 
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<n 


Prob (carrier  survives) 


•IT 


M 


If' 

ljf(0)  (i  -  .n<i-pA.  ) 


id 


N- 

i 


*TV 

i« 


(d> 

\f(d) 

UCj 


where. 


is  as  in  Eq'n  (2.5) , 


M 

(i  -  na-p  ) 

it i  3 1 


id 


3 


P^-  is  as  in  Eq’n  (2.6), 

(aj.bj)  are  the  limits  of  the  ith  AAW  sector,  and 
(c*,dj)  are  the  limits  of  the  j.th  ASW  sector. 


(5.1) 


5.2.2  DIFFERENT  PROBABILITY  DISTRIBUTIONS  OF  TARGETS 


The  introduction  of  different  probability  distributions  of  the  targets  on  a 
per  sector  basis,  would  add  some  flavour  to  the  model.  It  might  be  meaningful 
to  allow  the  target  threat  to  be  denser,  say,  in  the  middle  of  the  sector,  and 
wane  at  the  boundaries.  Then,  instead  of  having  a  single  f (d)  (a  uniform  dis¬ 
tribution)  for  every  sector,  we  would  have  a  series  of  distributions,  f,  (d) , 
f .(d),  . ,f,  (d) ,  where  G+H  =  total  number  of  (AAW  and  ASW)  sectors. 
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5.2.3  THE  CARRIER  AS  A  CONTRIBUTING  SHIP 


Our  model  precludes  the  assistance  of  the  carrier  itself.  In  a  realistic  situ¬ 
ation,  the  carrier  would  be  counted  on  to  provide  some  defense  of  its  own.  We 
could  furnish  the  carrier  with  its  own  AAW  (2^)  and  ASW  (2^)  capabilities.  The 
carrier  itself  would  then  be  treated  as  another  ship,  with  the  additional 
stipulation  that  its  location  not  be  changed;  i.e.  the  carrier  remains  at 
(0,0)  and  we  maximize  its  overall  survival  probability  only  over  the  positions 
of  the  AAW  and  the  ASW  ships,  as  before. 


5.2.4  INTRODUCTION  OF  GAME  THEORY 


Another  way  to  expand  this  problem  is  to  formulate  the  problem  in  terms  of  a 
game.  We  consider  two  commanders,  one  responsible  for  AAW  ships,  the  other  for 
ASW  ships,  in  the  spirit  of  the  Naval  CWC  command  -  and  -  control  doctrine. 


This  game-theoretic  approach  is  appealing.  Under  this  scheme,  each  commander 
has  a  local  objective  of  maximizing  the  probability  of  survival  of  the  carrier 
given  an  AAW  or  an  ASW  threat,  depending  on  his  warfare  responsibility.  They 
share  the  (common)  global  objective  of  maximizing  the  total  survival  probabil¬ 


ity  of  the  carrier,  given  both  types  of  threats. 
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The  commander's  decisions  of  where  to  place  his  ships  in  order  to  realize  his 
goal,  are  a  direct  result  of  maximizing  his  objective  function  with  respect  to 
the  positions  of  the  ships  which  are  in  his  control,  while  holding  the  posi¬ 
tions  of  the  other  commander's  ships  fixed. 

The  game  persists  as,  at  each  iteration,  each  commander  attempts  to  maximize 
his  objective  function,  given  the  previous  decision  of  the  other.  The  game 
terminates  when  we  converge  to  an  equilibrium.  This  is  Known  as  the  Nash 
strategy . 

Now,  we  state  the  conditions  which  guarantee  the  existence  of  a  Nash  solution. 

Necessary  and  Sufficient  Conditions  for  a  Nash  Solution: 

If  u,-6  U*»  U;  compact,  convex,  non-empty,  J;  (u(  ,ut, . ,un)  concave  in  Uj , 

continuous,  then  there  exists  a  Nash  equilibrium. 

We  formulate  these  ideas  mathematically.  Let  commander  l  be  the  commander 
responsible  for  the  AAW  threats,  and  hence  the  AAW  ships.  Let  J,  be  his  objec¬ 
tive  function.  Then, 
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J(  =  Prob (carrier  survives  given  an  AAW  threat) 


b; 

f(8) 

a, 


r 

(i  -  n(i-pA.  ))  d# 

I--4 


"»* 


(5.2) 


where , 

PA.  is  as  in  Eq'n  (2.5), 

(a;,bt-)  are  the  limits  of  the  ith  AAW  sector, 

N*  is  the  number  of  AAW  targets  in  the  ith  AAW  sector,  and 
G  is  the  number  of  AAW  sectors. 


Let  this  commander  have  M  ships,  with  coordinates:  (x, ,  y,  ) ,  (x^,yr) - - 

,  • ,  • 


Similarly,  let  commander  2  be  the  commander  responsible  for  the  ASW 
threats,  and  therefore  the  ASW  ships.  Let  be  his  objective  function. 

Then, 
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=  Proto (carrier  survives  given  an  ASW  threat) 


where, 

Pg.  is  as  in  Eq'n  (2.6), 

(c|,dj)  are  the  limits  of  the  ith  ASW  sector, 

Qi  is  the  number  of  ASW  targets  in  the  ith  ASW  sector,  and 
H  is  the  number  of  ASW  sectors. 


If  there  are  a  total  of  T  ships,  then  since  commander  1  has  M  ships,  commander 
2  has  (T  -  M)  ships,  with  coordinates:  (xil..t,y_.J  ,  (x  ,y  ),..,  , . .  (x^  ,y_). 

*  ""  M  *x  /Wrl_  T  ' 

Every  ship  has  both  AAW  and  ASW  capabilities,  and  hence  characteristic  lx.A-  , 
,  and  values.  Therefore,  any  ship,  whether  under  the  responsibil¬ 
ity  and  control  of  commander  l  or  commander  2,  contributes  to  both  J,  and  J4, 
just  as  it  did  in  the  centralized  case.  On  the  other  hand,  commander  l  only 
has  control  over  the  locations  of  his  (AAW)  ships.  Consequently,  he  maximizes 
J,  with  respect  to  the  positions  of  his  ships,  while  Keeping  the  positions  of 
the  ASW  ships  fixed.  Similarly  for  the  ASW  commander. 
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The  existing  computer  program  for  the  centralized  problem  can  be  used  as  a 
prototype  for  the  game.  The  latter  requires  the  interplay  of  two  such 
programs,  with  slight  modifications.  The  program  could  be  divided  into  two 
similar  modules,  one  for  each  commander.  The  module  for  commander  1  would  con¬ 
sist  of  the  following: 


1.  the  objective  function,  J( 

2.  an  array  containing  the  location  coordinates  of  his  ships 

3.  an  array  containing  the  location  coordinates  of  the  ASW  ships 

4.  a  facility  to  separate  the  PA|  ,  i  =  l,M  (the  a£w  ships)  from 
the  PS;  ,  i  =  M+1,T  (the  ASW  ships).  When  maximizing  J,  ,  it  is 
precisely  the  M  AAW  ships  and  their  effect  on  the  P^-  ,  i  =  l,M 
that  are  being  controlled.  The  (T  -  M)  ASW  ships  contribute  to 
J,  ,  but  their  positions  cannot  be  changed  by  the  AAW  commander 
These  positions  are  obtained  from  the  (previous)  decision  of  the 
ASW  commander  and  are  held  fixed  in  the  PA-  ,  i  =  M+l,T.  We  then 
write  J,  as: 
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rk' 

fA 

T 

] 

=TT 

i*-» 

f(S) 

J<X\ 

i 

(i  -  nu-p 

)  rm-pA.  ) 

|  ' 

)  dtf  J 

where. 


N; 


PA-  is  as  in  Eq'n  (2.5)  , 

N •  is  the  number  of  AAW  targets  in  the  ith  AAW  sector, 
(a*  ,  b- )  are  the  limits  of  the  ith  AAW  sector,  and 
G  is  the  number  of  AAW  sectors. 


(5.4) 


5.  an  optimization  and  integration  routine  as  in  the  original 
program 

6.  a  facility  to  pass  freely  the  solution  (i.e.  the  new  coordinates 
of  the  AAW  ships)  of  the  optimization  to  commander  2,  so  that  he 
may  proceed  with  the  maximization  of  . 

A  parallel  module  is  clearly  required  for  commander  2. 

It  is  convenient  to  present  the  game  in  terms  of  an  algorithm: 

Step  l.  Each  commander  provides  initial  conditions  for  the  positions  of 
his  ships. 
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Step  2.  Commander  1  maximizes  J,  with  respect  to  (x,  ,y(  )  ,  .... 

holding  ,yMM),  -  (xr  ,yr  )  fixed,  from  commander  2's 

previous  decision. 


Step  3.  Commander  2  maximizes  with  respect  to  (x^, y^)  »  • —  (x^,  ,yr^ 
holding  (x^y,),  ....  ( , y^)  fixed,  from  commander  l's  previous 
decision. 


Step  4.  IF 


new 

[(x,,y,),  . ..,  (xM,yH)j[  =  jjtx,  ,y,) ,  .  ..,(xM,yM)j 


new  c 

IF  {(x^y^) . (x,  ,yr  )]  =  {(x^y*^) . <xT  ,yr  \ 

THEN  STOP;  we  have  converged  to  a  Nash  strategy  and  an 


equilibrium  has  been  reached. 


The  optimal  values  of  J,  and  Jv  are  their  current 
values.  The  optimum  locations  of  the  ships  are  the 
current  values  of  (x,,yt),  . . . ,  (x^y^) ,  ( x  ,  yM  J  , . . , 
(xT  ,yr ) . 

ELSE  return  to  Step  2. 
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APPENDIX  A  PROOF  OF  CLAIM  3.1 


Proof  of  CLAIM  3.1: 

We  want  to  show  that  maximizing  the  survival  probability  of  the  carrier  with 
respect  to  the  locations  of  the  ships,  given  M  ships  and  1  target,  is  the  same 
as  maximizing  it  for  M  ships  and  N  targets.  We  consider  both  of  these  cases 
and  show  that  they  are  the  same. 


Case  1:  M  Ships,  1  Target 

We  wish  to  maximixe: 


max 

x;  '*.• 


(i 


/V\ 

-  nu- 

1-4 


i  =  1,..,M 


In  other  words,  we  want  to  solve  for  (x‘  ,y(-  ) ,  i  =  l,...,M 


-  J 

r  b 

7 

J«f(0) 

u  -  RU-  Pjfc- 

)) 

d0  J 

=  0 

a.  ( 

M 

(1  -  [1(1-  P*. 

f*f 

)) 

dfl  J 

=  0 

\Cb 

AH 

2 

j 

1  ifi9) 

U  -  -ilU-P*; 

)) 

d0  j 

=  0 

2_  1 

r  r*» 

M 

)) 

d0  J 

•  \  tie) 

A 

(l  -  FI  (  i-Pa; 

i*i 

=  0 
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*4 


^4$ 


f($)  (i  -  nd-pA.  ) 

Q  Pi  ' 


t(9)  (i  -  n(i-p^.  ) 

it.  71 1 


Case  2:  M  Ships,  N  Targets 


Here  we  want  to  maximize: 


tit)  (i  -  nu-pA.  ))  dtf 


x;  ,y; 


i  =  l ,  . .  .  ,  K 


In  other  words,  we  want  to  solve  for  (x,',y,),  i  = 


?  f  r-b  A 

a  -  nu-pA|.  )) 

(1-,C(1-pa;  }) 


(1  -  •Pil<1-pA.-  >>  66 

l;,Jfbf(0)  a  -  ,nu-pA-  ))  as 


is  f%<*>  a  -  J?U-Pa;  ))  Ad  >1 


Let  us  consider,  for  example: 

n  r  r  b 


i  tfi 


fid)  (i  -  nd-p*.  )) 

I  r-i 
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We  want  to  show  that  Eq'n  (A.l)  is  the  same  as  : 

3 


Jt  Ow 

ax 


a 


f (0)  d  -  n(i-p..  )) 

a  i-H  *' 


=  o 


(A. 2) 


SO, 


rrb 

M 

] 

Ua"*' 

(i  -  nu-pA.  )> 

i-i  M| 

ae  J 

=  0 


(A.l) 


"U 


dx; 


M 


fan  (i  -  jiu-pA.  )) 

i  t -«  M  < 


Tb 

\f(D)  d-.na-p.. 

Ja  "  * 


d0 


3 


))  <30 


j 


N-l 


(A. 3) 


In  order  to  show  that  Eq'n  (A. 3)  produces  the  same  solution  as  Eq'n  (A. 2) 
we  consider  3  possibilities. 


Either  i)  N  =  0 
or 

fb 

ii) If (0)  (1 

£rb 

iii)  V  (0)  (1 


or 


n< 

i -i 


nd-pA;  ))  d0 


nd-pAi.  ))  d0 

i=.  M ' 


=  0 


=  0 
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If  iii)  holds,  we  are  done.  So,  we  need  to  show  that  i)  and  ii)  cannot 
hold. 


i)  N  is  the  number  of  targets,  so  N  ^  o. 


ii)  If 


1 


M 


f(9)  (i-IKi-Pa-  ))  ad 

in 


*  0,  then 


either 

a) 

f  (9)  =  0 

or 

b) 

(a  ,b  )  has 

/* 

or 

c) 

(1  -  nu-p„; 

Then: 


a)  f(9)  is  the  uniform  distribution,  and  so  f(9)  /  0,  in  the  area 
of  interest,  i.e.  f{8 )  =  l/(b  -  a),  a  <=  9  $  b, 

b)  sector  (a  ,b  )  has  positive,  non-zero  area, 
and, 

* 

c)  (1  -  fl  (1-Pa-  ))  *  0  since  P.,+1  VI  . 

m  A, 


Therefore,  iii)  holds  always,  and  Case  l  and  Case  2  produce  the  same 
solution,  and  the  optimal  ship-to-carrier  distance,  R,  is  independent 
of  the  number  of  (independent  and  identically  distributed)  targets. 
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APPENDIX  B  ABOUT  THE  PROGRAM 


B.l  HOW  THE  MODEL  IS  IMPLEMENTED 


This  section  deals  with  the  numerical  implementation  of  the  model.  It 


explains  how  the  computer  program  works. 


All  experiments  were  run  at  the  Massachusetts  Institute  of  Technology  on  an 


I . B.M.  370/168. 


The  cost  of  a  typical  optimization  run  in  CPU  time  is  6  seconds.  The  associ¬ 
ated  dollar  cost  is  $  1.25. 


B.l.l  THE  MAIN  PROGRAM,  LIBRARY  ROUTINES  AND  SUBROUTINES 


The  program  makes  use  of  several  modules  and  subroutines.  The  main  program  has 
three  tasks.  It 

1)  calls  an  optimization  routine  (SUMSNO)  to  maximize  the 
objective  function, 

2)  calls  a  subroutine  (READ3)  to  read  the  input  data,  and 


3)  prints  the  final  positions  of  the  ships. 
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SUMSNO  (5)  is  a  NONLIN  library  subroutine  which  minimizes  general  uncon¬ 
strained  objective  functions  of  low  enough  order.  We  use  it  to  find  the  opti¬ 
mal  locations  of  the  ships  -  i.e.  those  coordinates  of  the  ships  that  maximize 
the  probability  of  survival  of  the  carrier.  At  each  iteration,  SUMSNO  "tries 
out"  a  set  of  ships'  coordinates  and  evaluates  the  objective  function,  i.e., 
the  survival  probability  of  the  carrier.  It  uses  an  approximation  to  both  the 
gradient  and  the  Hessian. 


We  recall  from  Chapter  2,  Section  2. 3. 2.1  Case  4,  that  the  most  general 
expression  for  the  survival  probability  of  the  carrier  in  a  sector  (a,b)  is: 


(B.l) 


where , 

M  =  number  of  ships, 

N  =  number  of  AAW  targets,  and 
Q  =  number  of  ASW  targets 

fid)  is  the  uniform  distribution  of  a  target  over  a  sector  (a,b) 
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P^.  is  the  Prob(ship  i  kills  an  AAW  target) 

Pjt-  is  the  Prot>(ship  i  kills  an  ASW  target) 

Since  the  probability  in  Eq'n  (B.l)  is  the  product  of  integrals,  we  need  a 
subroutine  that  will  evaluate  the  integrals  numerically.  D01BAF,  a  subroutine 
which  is  part  of  the  NAG  library,  is  used  to  do  this.  We  invoke  D01BAF,  the 
integration  subroutine,  at  each  iteration,  from  QDRTF ,  the  subroutine  that 
evaluates  the  objective  function.  There  are  two  function  subprograms,  AAWFNl 
and  ASWFNl  (for  AAW  and  ASW  sectors  respectively)  which  are  called  by  D01BAF 
each  time  it  is  necessary  to  integrate  over  a  particular  sector. 

Please  refer  to  the  References  (5)  and  (6)  for  more  information  on  SUMSNO  and 
DOlBAF  respectively. 

B.l. 2  THE  DATA  FILE 

The  data  is  stored  in  a  data  file  called  KSH  DT.  One  modifies  the  data  by 
editing  the  contents  of  KSH  DT.  The  subroutine  that  reads  the  input,  READ3, 
assumes  a  fixed  line-by-line  format  for  the  input  file.  This  implies  that 
maintaining  the  format  is  crucial  in  order  for  the  program  to  work  correctly. 
A  sample  format  of  the  data  file  is  included. 
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This  report  lends  itself  to  experiments  in  which  the  number  of  AAW  (ASW)  ships 
and/or  sectors  and/or  targets  is  varied.  The  data  sits  in  arrays  whose  dimen¬ 
sions  may  be  easily  altered.  This  requires  adding  (for  augmenting  the  size  of 
an  array)  or  deleting  (for  reducing  the  size  of  an  array)  lines  from  KSH  DT. 
It  is  in  this  situation  that  one  must  pay  attention  to  the  format.  Please 
refer  to  the  comments  in  the  subroutine  READ3  of  the  program  for  any  further 
information.  The  computer  program  can  be  found  in  Appendix  C. 


B.1.3  THE  ITERATION  SUMMARY 


SUMSNO  supplies  an  iteration  summary  via  an  internal  subroutine  ITSUMU. 
ITSUMU  provides  basic  information  about  the  objective  function,  such  as  its 
current  value  at  a  particular  iteration.  It  provides  the  optimal  solution  in 
Cartesian  coordinates.  Our  computer  program  also  gives  the  optimal  solution  in 
polar  coordinates. 

Since  the  commands  to  run  the  program  are  coded  in  an  EXEC  file,  it  is  easy  to 
use  the  program.  The  EXEC  command  is: 

DFKG3 
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For  a  more  detailed  description  of  the  program  and  an  explanation  of  its 
parameters,  we  advise  the  reader  to  consult  the  program  itself.  It  is  docu¬ 
mented  extensively  and  may  be  found  in  Appendix  C. 
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APPENDIX  C  THE  COMPUTER  FROGRAM  AND  ITS  DOCUMENTATION 


C.l  THE  COMPUTER  PROGRAM  AND  ITS  DOCUMENTATION 


The  actual  program  follows. 
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*  *  *  PURPOSE  OF  THIS  PROGRAM  *  *  * 

THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  MAXIMIZE  THE  PROBABILITY  OF  SURVIVAL 
OF  A  CARRIER,  GIVEN  AIR  AND  SUB  THREATS.  THE  CARRIER  HAS  NO  DEFENSIVE 
CAPABILITIES  OF  ITS  OWN.  IT  IS  DEFENDED  BY  AAW  (ANTI-AIR  WARFARE)  AND 
ASW  (ANTI-SUB  WARFARE)  SHIPS.  THE  USER  SUPPLIES  INPUT  CONCERNING 
THE  NUMBER  OF  SHIPS,  THE  CHARACTERISTIC  PARAMETERS  OF  THE  SHIPS,  THE 
NUMBER,  TYPE  (AAW  OR  ASW)  AND  LOCATION  OF  THE  THREAT  SECTORS,  AND  THE 
NUMBER  OF  TARGETS.  ALL  THIS  IS  EXPLAINED  IN  GREATER  DETAIL  AS  IT  IS 
ENCOUNTERED  IN  THE  PROGRAM. 

«  «  «  TO  RUN  THE  PROGRAM  *  *  * 

IN  ORDER  TO  RUN  THIS  PROGRAM,  ONE  USES  THE  EXEC  FILE  CALLED: 

1)  DFKG3 ,  FOR  PRINTINGS  OF  INITIAL  INPUT  AND  FINAL  OUTPUT  ONLY, 

IT  IS  INVOKED  SIMPLY  BY  TYPING:  DFKG3 

*  *  «  GENERAL  OVERVIEW  «  *  » 

THE  PROGRAM  IS  DIVIDED  INTO  FIVE  MODULES.  THEY  ARE: 

1.  MAIN  PROGRAM: 

THE  MAIN  PROGRAM  HAS  FIVE  FUNCTIONS.  IT 

I.  CALLS  READ3 ,  THE  SUBROUTINE  THAT  READS  IN  A  THE  INPUT  PARAMETERS. 
II.  CONVERTS  DEGREES  TO  RADIANS 
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.  PRINTS  INFORMATION  CONCERNING  THE  INPUT  PARAMETERS 

.  CALLS  SUMSNO ,  THE  OPTIMIZATION  ROUTINE 

.  COMPUTES  (VIA  SUMSNO)  THE  OPTIMAL  SHIP-TO-CARRIER  DISTANCES,  R ( I ) 
OF  EACH  SHIP  I,  AND  THEIR  ASSOCIATED  OPTIMAL  ANGLES.  IT  THEN 
PRINTS  THEM  AT  THE  END  OF  THE  OUTPUT. 

SUBROUTINE  QDRTF 

.  EVALUATES  THE  OBJECTIVE  FUNCTION,  F;  F  =  SURVIVAL  PROBABILITY  OF 
THE  CARRIER.  SINCE  SUMSNO  ACTUALLY  MINIMIZES  F,  AND  WE  WANT  TO 
MAXIMIZE  OUR  OBJECTIVE  FUNCTION,  WE  MINIMIZE  (-  F) . 

.  SINCE  F  IS  THE  PRODUCT  OF  INTEGRALS,  WE  NEED  AN  INTEGRATION 
ROUTINE  TO  EVALUATE  THE  INTEGRAL (S) .  THIS  IS  DONE  BY  D01BAF. 
DOIBAF  IS  CALLED  BY  QDRTF. 

.  CALLS  A  FUNCTION  SUBROUTINE,  AAWFN1 ,  TO  EVALUATE  THAT  PART 
OF  F  ASSOCIATED  WITH  AAW  THREATS. 

.  QDRTF  CALLS  A  FUNCTION  SUBROUTINE,  ASWFN1,  TO  EVALUATE  THAT  PART 
OF  F  ASSOCIATED  WITH  ASW  THREATS.  SINCE  AAW  THREATS  ARE 
INDEPENDENT  OF  ASW  THREATS,  F  =  -  RR  »  SS, 

WHERE, 

RR  =  PROB (CARRIER  SURVIVES  GIVEN  AAW  THREATS),  AND 
SS  =  PROB (CARRIER  SURVIVES  GIVEN  ASW  THREATS) 

SUBROUTINE  AAWFN1 : 

IT  EVALUATES  THAT  PART  OF  THE  SURVIVAL  PROBABILITY  OF  THE  CARRIER 
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C  DUE  TO  AAW  THREATS. 

C 

C 

C  4.  SUBROUTINE  ASWFN1: 

C 

C  IT  EVALUATES  THAT  PART  OF  THE  SURVIVAL  PROBABILTI Y  OF  THE  CARRIER 
C 

C  DUE  TO  ASW  THREATS. 

C 

C 

C  5.  SUBROUTINE  READ 3 : 

C 

C  IT  READS  IN  ALL  THE  INPUT  PARAMETERS. 

C 

C 

C 

C 

C 

C  THE  FOLLOWING  ARE  THE  PARAMETERS  THAT  APPEAR  IN  THIS  PROGRAM. 

C 

C  THEY  APPEAR  IN  ALPHABETICAL  ORDER. 

C 

C 

£XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 

C  PARAMETER  USAGE 

C 

QXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

c 

C  X  X  x  INPUT  PARAMETERS  *  *  * 


2 


9 


C  AAWEXP ( I ) -  THE  NUMBER  OF  TARGETS  IN  THE  ITH  AAW  THREAT  SECTOR. 

C 

C  AAWLL(I) .  THE  LOWER  LIMIT  OF  THE  ITH  AAW  THREAT  SECTOR. 

C 

C  AAWUL(I) .  THE  UPPER  LIMIT  OF  THE  ITH  AAW  THREAT  SECTOR. 

C 

C  ASIGMA(I) _  THE  AAW  CAPABILITY  FACTOR  OF  THE  ITH  SHIP. 

C 

C  ASWEXP ( I ) -  THE  NUMBER  OF  TARGETS  IN  THE  ITH  ASW  THREAT  SECTOR. 

C 

C  ASWLL(I) .  THE  LOWER  LIMIT  OF  THE  ITH  ASW  THREAT  SECTOR. 

C 

C  ASWUL(I) .  THE  UPPER  LIMIT  OF  THE  ITH  ASW  THREAT  SECTOR. 


Appendix  C  The  Computer  Program  and  its  Documentation 


I 

PAGE 


C 

C  G .  THE  NUMBER  OF  AAW  THREAT  SECTORS. 

C 

C  H .  THE  NUMBER  OF  ASW  THREAT  SECTORS. 

C 

C  KA ( I ) .  AN  AAW  PARAMETER  REFLECTING  THE  EFFECT  OF  THE  ITH  SHIP'S 

C 

C  DISTANCE  TO  THE  CARRIER. 

C 

C  KS ( I ) .  AN  ASW  PARAMETER  REFLECTING  THE  EFFECT  OF  THE  ITH  SHIP'S 

C 

C  DISTANCE  TO  THE  CARRIER. 

C 

C  M .  THE  TOTAL  NUMBER  OF  SHIPS  (AAW  AND  ASW). 

r* 

W 

C  NN .  THE  TOTAL  NUMBER  OF  SHIP  COORDINATES  (TWICE  THE  NUMBER 

C 

C  OF  SHIPS) . 

C 


C  SSIGMA ( I ) -  THE  ASW  CAPABILITY  FACTOR  OF  THE  ITH  SHIP. 

C 

C 

C 

C 

c 

C  *  *  *  INTERNAL  ARRAYS  REQUIRED  BY  SUMSNO  *  *  * 

C 

C 

C  WE  GIVE  A  BRIEF  DESCRIPTION  OF  THE  ARRAYS  WHICH  SUMSNO  USES 
C 

C  INTERNALLY.  FOR  A  MORE  DETAILED  DESCRIPTION,  AND  FOR  MORE  INFORMATION 
C 

C  ABOUT  SUMSNO  AND  ITS  ARRAYS  AND  SUBROUTINES,  PLEASE  REFER  TO  THE 
C 

C  REFERENCES  OF  THIS  REPORT. 


C 

C 

C  D .  (INPUT/OUTPUT)  A  SCALE  VECTOR  SUCH  THAT  D(I)*X(I), 

C 

C  1  =  1, 2 - - N  ARE  ALL  IN  THE  COMPARABLE  UNITS.  SINCE 

C 

C  WE  REQURE  NO  SCALING,  D  IS  INITIALIZED  TO  ONES.  THE 

r* 

C  DEFAULTS  PROVIDED  BY  DFALTU  (SEE  BELOW)  REQUIRE  THE 


CALLER  TO  SUPPLY  D. 


176 
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C  DFALTU (IV , V) .  (INPUT)  A  SUBROUTINE  OF  INPUT  DEFAULT  VALUES.  SEE  BELOW. 
C 

C  IV .  (INPUT/OUTPUT)  AN  INTEGER  VALUE  ARRAY  OF  LENGTH  AT 

C 

C  LEAST  39  THAT  HELPS  CONTROL  THE  SUMSNO  ALGORITHM.  IT  IS 

C 

C  ALSO  USED  TO  STORE  SEVERAL  INTERMEDIATE  QUANTITIES. 

C 

C  N .  (INPUT)  THE  NUMBER  OF  VARIABLES  ON  WHICH  F  DEPENDS, 

C 

C  I.E.  THE  NUMBER  OF  COMPONENTS  IN  X. 

C 

C  QDRTF .  (INPUT)  A  SUBRROUTINE  THAT,  GIVEN  X,  CALCULATES  F(X). 

C 

C  QDRTF  MUST  BE  DECLARED  EXTERNAL  IN  THE  CALLING  PROGRAM. 

C 

C  IT  IS  INVOKED  BY 

C 

C  CALL  QDRTF (N ,  X,  NF,  F,  UIP,  URP,  UFP) 

C 

C  NF  IS  THE  INVOCATION  COUNT  FOR  QDRTF.  THE  OTHER 

C 

C  PARAMETERS  ARE  AS  DESCRIBED  ABOVE  AND  BELOW. 

C 

C  SUMSNO .  MINIMIZE  (MAXIMIZE)  GENERAL  UNCONSTRAINED  OBJECTIVE 

C 

C  FUNCTION  USING  FINITE-DIFFERENCE  GRADIENTS  AND  SECANT 

C 

C  HESSIAN  APPROXIMATIONS. 

C 

C  X .  (INPUT/OUTPUT)  BEFORE  (INITIALLY)  CALLING  SUMSNO,  THE 

C 

C  CALLER  SHPULD  SET  X  TO  AN  INITIAL  GUESS  AT  X*.  WHEN 

C 

C  SUMSNO  RETURNS,  X  CONTAINS  THE  BEST  POINT  SO  FAR 

C 

C  FOUND . 

C 

C  UFP .  (INPUT)  USER  EXTERNAL  SUBROUTINE  OR  FUNCTION  PASSED 

C 

C  WITHOUT  CHANGE  TO  QDRTF. 

C 

C  UIP .  (INPUT)  USER  INTEGER  PARAMETER  ARRAY  PASSED  WITHOUT 

C 

C  CHANGE  TO  QDRTF. 
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C 

C  URP . 

C 

c 

c 

C  V . 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

C  ANGLE (I) 
C 

C  RAD ( I ) . . 
C 
C 
C 

c 

c 

c 

c 

c 

C  C,  E,  I, 
C 

C  DATAN... 
C 

C  DSQRT... 
C 

C  MINNUM.. 
C 
C 
C 

C  PHI (I) .. 
C 
C 
C 

C  PI . 


(INPUT)  USER  FLOATING-POINT  PARAMETER  ARRAY  PASSED 
WITHOUT  CHANGE  TO  QDRTF. 

(INPUT/OUTPUT)  A  FLOATING-POINT  VALUE  ARRAY  OF  LENGTH 
AT  LEAST  67  +  (N«(N+l9)/2)  THAT  HELPS  CONTROL  THE 
SUMS NO  ALGORITHM.  IT  IS  ALSO  USED  TO  STORE  SEVERAL 
INTERMEDIATE  QUANTITIES. 


*  *  *  OUTPUT  PARAMETERS  *  *  * 

THE  OPTIMAL  ANGLE  OF  THE  ITH  SHIP. 

THE  OPTIMAL  SHIP-TO-CARRIER  DISTANCE  OF  THE  ITH  SHIP. 


*  *  *  MISCELLANEOUS  *  *  « 

INTEGERS  USED  FOR  INDEXING  PURPOSES. 

DOUBLE  PRECISION  ARCTANGENT  FUNCTION. 

DOUBLE  PRECISION  SQUARE  ROOT  FUNCTION. 

IMPLEMENTATION  DEPENDENT  CONSTANT.  THE  SMALLEST 
EXPONENT. 

THE  OPTIMAL  ANGLE  OF  THE  ITH  SHIP.  IT  IS  USED  IN  THE 
PROGRAM,  BUT  NOT  PRINTED  OUT. 

THE  IRRATIONAL  NUMBER,  3.141592653589793 
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C 

C  SMALL .  IMPLEMENTATION  DEPENDENT  CONSTANT.  THE  SMALLEST 

C 

C  NUMBER  THAT  CAN  BE  EVALUATED. 

C 

C  Y  (I) .  ITS  THE  VALUE  OF  Y(I)/X(I),  WHERE,  (X(I),Y(I))  ARE  THE 

C 

C  OPTIMAL  CARTESIAN  COORDINATES  OF  THE  ITH  SHIP.  IT  IS 

C 

C  USED  IN  THE  PROGRAM,  BUT  IS  NOT  PRINTED  OUT. 


C 

C 

C 

c 

C  *  *  *  DECLARATION  OF  VARIABLES  *  *  * 

C 

C 

EXTERNAL  DFALTU ,  QDRTF ,  SUMSNO 
INTEGER  I,  IV (39) ,  UIP(l) 

DOUBLE  PRECISION  D(20) ,  URP(20,3),  V(457),  X(20) 

DOUBLE  PRECISION  KA (150) , KS (150) ,  ASIGMA(150),  SSIGMA(150) 
DOUBLE  PRECISION  AAWEXP(150),  ASWEXP(150) 

DOUBLE  PRECISION  PHI(150),  RAD(150) 

DOUBLE  PRECISION  DATAN,  DSQRT 
DOUBLE  PRECISION  Y(150),  ANGLE (150) 

DOUBLE  PRECISION  AAWLL(150),  AAWUL(150),  ASWLL(lSO) ,  ASWUL(150) 
REAL* 8  PI 
REAL *8  MINNUM 
REAL *8  ''MALL 

INTEGER  C,  E,  G,  H,  M,  NN 
COMMON  /V/  KA 
COMMON  /W/  KS 
COMMON  /U/  ASIGMA 
COMMON  /W/  SSIGMA 
'  COMMON  /C/  AAWEXP 
COMMON  /D/  ASWEXP 
COMMON  /Y/  AAWLL,  AAWUL 
COMMON  /Z/  ASWLL,  ASWUL 
COMMON  /ASECT/  G 
COMMON  /SSECT/  H 
COMMON  /SHIP/  M 
COMMON  /X/  X 
COMMON  /MIN/  MINNUM 
COMMON  /SM/  SMALL 
PI  =  3.141592653589793 
C 
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*  *  *  HOW  THE  VALUES  ARE  READ  IN  THE  DATA  FILE  *  *  * 


THE  VALUES  ARE  READ  IN  THE  FOLLOWING  ORDER: 


MINNUM 
NN 
D  ( I ) 

X(I) 

M 

KA(I) 
KS(I) 
ASIGKA (I) 
SSIGMA ( I ) 
G 

AAWLL(I) 

AAWUL(I) 

AAWEXP(I) 

H 

ASWLL(I) 
ASWUL(I) 
ASWEXP ( I ) 


SMALL 


*  *  *  CALL  READ 3  TO  GET  THE  INPUT  FROM  THE  DATA  FILE  «  *  * 


CALL  READ3(G,H,M,NN,D,X,KA,KS,ASIGMA, SSIGMA, AAWLL , AAWUL , AAWEXP , 
1ASWLL , ASWUL , ASWEXP , MINNUM , SMALL ) 


*  *  »  CONVERT  SECTOR  DEGREES  TO  RADIANS  *  *  « 

THE  LIMITS  OF  THE  SECTORS  WERE  GIVEN  IN  DEGREES.  WE  NEED  TO  GIVE 
THE  LIMITS  IN  RADIANS.  THE  NEXT  TWO  "DO  LOOPS"  CONVERT  THE  DEGREES 
TO  RADIANS. 
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C 

c 

c 


51 

c 


52 


DO  51  I  =  1,G 

AAWLL(I)  =  AAWLL(I)  *  PI/180. D+0 
AAWUL(I)  =  AAWUL(I)  *  PI/180. D+0 
CONTINUE 


DO  52  I  =  1,H 

ASWLL(I)  =  ASWLL ( I )  »  PI/180. D+0 
ASWUL(I)  =  ASWUL(I)  *  PI/180. D+0 
CONTINUE 


C 

C 

C 

C 

C 

C 

C 


400 


16 


*  *  *  PRINT  VARIOUS  INFORMATION  «  *  * 

WRITE  (6,400) 

FORMAT  (IX,  'KA(I)  IS  AAW  PARAMETER  OF  ITH  SHIP  REFLECTING’, 
1*  EFFECT  OF  ITS  DISTANCE  TO  CARRIER') 

WRITE  (6,16)  (I,  KA (I ) ,  I  =  1,M) 

FORMAT  (IX,  'KA( ' , 12, ' )  *  ',  F10.5) 


WRITE  (6,401) 

401  FORMAT  ('O',  'KS(I)  IS  ASW  PARAMETER  OF  ITH  SHIP  REFLECTING', 
1'  ',' EFFECT  OF  ITS  DISTANCE  TO  CARRIER') 

WRITE  (6,161)  (I,  KS ( I ) ,  I  =  1,M) 

161  FORMAT  (IX,  'KS(',I2,')  =  ',  F10.5) 


WRITE  (6,402) 

402  FORMAT  ('O',  'ASIGMA(l)  IS  THE  AAW  CAPABILITY  FACTOR  OF  SHIP  I') 
WRITE  (6,17)  (I,  ASIGMA ( I ) ,  I  =  1,M) 

17  FORMAT  (IX,  ' ASIGMA ( ' ,  12,')  =  ’,  F10.5) 


WRITE  (6,403) 

403  FORMAT  ('O',  'SSIGMA(I)  IS  THE  ASW  CAPABILITY  FACTOR  OF  SHIP  I') 
WRITE  (6,18)  (I,  SSIGMA ( I ) ,  I  =  1,M) 

18  FORMAT  (IX,  'SSIGMA(',  12,’)  =  ',  F10.5) 


WRITE  (6,404) 
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404  FORMAT  ('O’,  'AAWLL(I)  IS  THE  LOWER  LIMIT  OF  THE  ITH  AAW  SECTOR') 
WRITE  (6,19)  (I,  AAWLL(I) ,  I  *  1,G) 

19  FORMAT  (IX,  ' AAWLL ( ' ,12, 1 )  =  ' ,  F10.5) 


WRITE  (6,405) 

405  FORMAT  ('O’,  ’AAWUL(I)  IS  THE  UPPER  LIMIT  OF  THE  ITH  AAW  SECTOR') 
WRITE  (6,20)  (I,  AAWUL(I),  I  =  1,G) 

20  FORMAT  (IX,  ' AAWUL ( ' , 12 , ' )  =  ',  F10.5) 


WRITE  (6,406) 

406  FORMAT  ('O',  'ASWLL(I)  IS  THE  LOWER  LIMIT  OF  THE  ITH  ASW  SECTOR') 
WRITE  (6,21)  (I,  ASWLL(I),  I  =  1,H) 

21  FORMAT  (IX,  ’ ASWLL ( ' , 12 , ’ )  =  ’,  F10.5) 


WRITE  (6,407) 

407  FORMAT  ('O',  'ASWUL(I)  IS  THE  UPPER  LIMIT  OF  THE  ITH  ASW  SECTOR’) 
WRITE  (6,22)  (I,  ASWUL(I),  I  =  1,H) 

22  FORMAT  (IX,  ' ASWUL ( ' , 12 , ' )  =  ',  F10.5) 


WRITE  (6,408) 

408  FORMAT  ('O',  'AAWEXP(I)  IS  THE  NUMBER  OF  TARGETS  IN  AAW  SECTOR  I') 
WRITE  (6,23)  (I,  AAWEXP ( I ) ,  I  =  1,G) 

23  FORMAT  (IX,  ' AAWEXP (' 12 ,' )  =  ’,  F10.5) 


WRITE  (6,409) 

409  FORMAT  ('O',  ’ASWEXP(I)  IS  THE  NUMBER  OF  TARGETS  IN  ASW  SECTOR  I') 
WRITE  (6,24)  (I,  ASWEXP ( I ) ,  I  =  1,H) 

24  FORMAT  (IX,  ' ASWEXP (' 12 ,' )  =  ',  F10.5) 


WRITE  (6,25) 

25  FORMAT  (16HONOPNGN  ON  QDRTF) 
IV(1)=0 


WE  USE  THE  SUMSNO  SUBROUTINE  WHICH  MINIMIZES  (MAXIMIZES)  A  GENERAL 
UNCONSTRAINED  OBJECTIVE  FUNCTION  USING  A  GRADIENT 
AND  A  HESSIAN  APPROXIMATION.  THE  MAIN  PROGRAM  REPEATEDLY 
MINIMIZES  (MAXIMIZES)  THE  OBJECTIVE  FUNCTION. 
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I 

C  THE  ROUTINE  QDRTF  EVALUATES  THE  OBJECTIVE  FUNCTION  F(X) 

C  DESCRIBED  IN  THE  MAIN  PROGRAM. 

C 
C 
C 

C  *  *  *  CALL  SUMSNO  TO  MAXIMIZE  THE  OBJECTIVE  FUNCTION  *  *  * 

C 
C 

CALL  SUMSNO (NN ,  D,  X,  QDRTF,  IV,  V,  UIP,  URP,  QDRTF) 

C 
C 
C 

Qx*«*****x****xx*xx*xxft*xx**x*x****x*****ft**ft*ft**x*x****x***x**xxx***xx* 

c 

C  WE  NOW  COMPUTE  THE  ANGLES  OF  THE  LOCATIONS  OF  THE  SHIPS. 

C  "PHI (I)"  IS  THE  ANGLE  ASSOCIATED  WITH  SHIP  I  AT  ITS  FINAL 

C  (OPTIMAL)  LOCATION. 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

c  WE  COMPUTE  THE  ANGLE  AS  FOLLOWS:  LET  SHIP  I  HAVE  FINAL  COORDINATES 
C  (X (E) ,X (F) ) .  THEN  THE  ANGLE  PHI (I)  ASSOCIATED  WITH  SHIP  I'S 

C  FINAL  (OPTIMAL)  POSITION  IS  GIVEN  BY  PHI (I)  =  ARCTAN (X (F) /X (E) ) , 

C  WHICH  IN  THIS  PROGRAM  IS  WRITTEN  AS  PHI(I)  =  DATAN ( Y ( I ) ) .  THIS 
C  COMPUTATION  RETURNS  PHI (I)  IN  RADIANS.  WE  WOULD  LIKE  TO  HAVE 
C  PHI (I)  IN  DEGREES,  FOR  SIMPLICITY  SAKE.  IN  ORDER  TO  CONVERT  TO 
C  DEGREES,  WE  NEED  TO  CONSIDER  ALL  POSSIBLE  COMBINATIONS  OF 

C  (X ( E) ,X(F) )  -  I.E.  ALL  QUADRANTS.  THIS  IS  SO  BECAUSE,  FOR  EXAMPLE, 

C  ARCTAN  (I.E.  DATAN)  MAKES  NO  DISTINCTION  BETWEEN  (X(F)/X(E))  AND 
C  (-X(E)/-X(F) ) ,  X(E) ,  X (F)  >  0.  SO  IT  DOES  NOT  DIFFERENTIATE 
C  BETWEEN  QUADRANTS  1  AND  3,  AND  QUADRANTS  2  AND  4. 

C  THIS  EXPLAINS  ALL  THE  "IF"  STATEMENTS. 

C 
C 

PI  =  3.141592653589793 


DO  94  I  =  1 , M 

C  =  (2  *  I)  -  1 
E  =  (2  *  I) 

Y ( I )  =  X(E)/X(C) 

PHI (I)  =  DATAN ( Y ( I ) ) 

IF (X (C)  .GE.  O.D+O  .AND.  X(E)  .GE.  O.D+O)  ANGLE ( I )  =  (180. D+0 
1*  PHI ( I) ) /PI 

IF (X (C)  .LE.  O.D+O  .AND.  X(E)  .GE.  O.D+O)  ANGLE ( I )  =  ((180. D+0 
1*  PHI  (I ) ) /PI )  +  180. D+0 

IF(X(C)  .GE.  O.D+O  .AND.  X(E)  .LE.  O.D+O)  ANGLE ( I )  =  ((180. D+0 
2*  PHI ( I ) ) /PI )  +  360. D+0 
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IF (X(C)  .LE.  O.D+O  .AND.  X(E)  .LE.  O.D+O)  ANGLE ( I )  =  ((180.D+0 
3*  PHI ( I ) ) /PI )  +  180. D+0 
WRITE (6, 90)  I,  ANGLE (I) 

90  FORMAT (IX,  ' ANGLE ( ' ,  12,  ')  =  F9.3,  3X,  ’DEGREES') 

94  CONTINUE 
C 
C 

C  WE  COMPUTE  THE  RADIUS  OF  EACH  OF  THE  M  SHIPS,  RAD (I). 

C 

DO  103  I  =  1 ,M 

C  =  (2  *  I)  -1 
E  ■  (2  *  I) 

RAD ( I )  =  DSQRT  ((X(C)  *  X(C))  +  (X(E)  *  X(E) ) ) 

WRITE  (6,30)  I,  RAD ( I ) 

30  FORMAT  (IX,  ' RAD ( ' ,  12,')  =  ',  F10.5) 

103  CONTINUE 
C 

STOP 

END 

C 

C 

C 

C 

C 

QXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C  QDRTF 

c 

Q*  ************  *  A  ********************  X  *******************  X  **************  * 

c 

C  SUBROUTINE  QDRTF  IS  THE  ONE  THAT  EVALUATES  THE  OBJECTIVE  FUNCTION 

C  F (X)  . 

C  SINCE  F (X)  IS  AN  INTEGRAL  ,  WE  NEED  A  SUBROUTINE  TO  ESTIMATE  THE 

C  VALUE  OF  THE  INTEGRAL.  THIS  IS  DONE  BY  D01BAF. 

C  MOST  OF  THE  PARAMETERS  WERE  EXPLAINED  AT  THE  BEGINNING  OF  THE 

C  PROGRAM.  ANY  NEW  PARAMETERS  APPEARING  HERE  ARE  EXPLAINED  BELOW. 

C 

Q*  ******  It  *  *  *  it  *  *  it  *****  it  *******  it  ****************************************  * 

c 

C  PARAMETER  USAGE 

C 

C*X*ft*ft******************ft*ft*ft**********ftftftftft****ft**X****»*****ft******* 

c 

c 

C  AAWANS ( I ) . . . .  THE  VALUE  OF  THE  INTEGRAL  DUE  TO  A  SINGLE  THREAT  IN  THE 
C 
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C  ITH  AAW  THREAT  SECTOR,  IS  RETURNED  AS  AAWANS(I). 

C 

C  AAWFN1 .  A  REAL  FUNCTION  CALLED  BY  D01BAF  TO  EVALUATE  THE 

C 

C  INTEGRAND. 

C 

C  AS  A  RESULT,  IT  IS  ALSO  THE  INTEGRAND;  I.E  THE  SURVIVAL 

C 

C  PROBABILITY  OF  THE  CARRIER  DUE  TO  A  THREAT  IN  THE  ITH 

C 

C  AAW  SECTOR. 

C 

C  AAWTOT ( I ) . . . .  THE  SURVIVAL  PROBABILITY  OF  THE  CARRIER  DUE  TO  ANY 
C 

C  NUMBER  OF  THREATS  IN  THE  ITH  AAW  SECTOR.  IT  IS 

C 

C  DETERMINED  AS: 

C 

C  AAWTOT ( I )  =  AAWANS(I)  **  AAWEXP(I), 

C 

C  WHERE,  AAWEXP(I)  IS  THE  NUMBER  OF  TARGETS  IN  THE  ITH 

C 

C  AAW  SECTOR. 

C 

C  ASWANS(I) . . . .  THE  VALUE  OF  THE  INTEGRAL  DUE  TO  A  SINGLE  THREAT  IN  THE 
C 

C  ITH  ASW  SECTOR,  IS  RETURNED  AS  ASWANS(I), 

C 

C  ASWFN1 .  A  REAL  FUNCTION  CALLED  BY  DOIBAF  TO  EVALUATE  THE 

C 

C  INTEGRAND. 

C 

C  AS  A  RESULT,  IT  IS  ALSO  THE  INTEGRAND;  I.E.  THE 

C 

C  SURVIVAL  PROBABILITY  OF  THE  CARRIER  DUE  TO  A  THREAT  IN 

C 

C  THE  ITH  ASW  SECTOR. 

C 

C  ASWTOT(I) . . . .  THE  SURVIVAL  PROBABILTIY  OF  THE  CARRIER  DUE  TO  ANY 
C 

C  NUMBER  OF  THREATS  IN  THE  ITH  ASW  SECTOR.  THE  VALUE 

C 

C  IS  DETERMINED  AS: 

C 

C  ASWTOT ( I )  =  ASWANS (I)  **  ASWEXP(I), 

C 
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C  WHERE,  ASWEXPC)  IS  THE  NUMBER  Sr  TARGETS  IN  THE  ITH 

C 

C  ASW  SECTOR. 

C 

C  DEXP .  DOUBLE  PRECISION  EXPONENTIAL  FUNCTION. 

C 

C  DOIBAF .  A  FUNCTION  THAT  COMPUTES  AN  ESTIMATE  OF  THE  DEFINITE 

C 

C  INTEGRALS,  AAWFN1  AND  ASWFN1. 

C 

C  D01BA2 .  A  SUBROUTINE  PROVIDED  BY  THE  NAG  LIBRARY  USED  TO 

C 

C  EVALUATE  AN  INTEGRAL  ON  A  FINITE  INTERVAL. 

C 

C  IFAIL .  THIS  IS  AN  ERROR  INDICATOR. 

C 

C  II .  USED  FOR  INDEXING  PURPOSES. 

C 

C  NF .  THE  NUMBER  OF  FUNCTION  EVALUATIONS. 

C 

C  NSTOR(l) .  THE  NUMBER  OF  ABSCISSAE  USED  IN  THE  EVALUATION  OF  THE 

C 

C  INTEGRAL. 

C 

C  RR .  VALUE  OF  THE  SURVIVAL  PROBABILITY  OF  THE  CARRIER  DUE 

C 

C  TO  AAW  THREATS  FROM  ALL  THE  AAW  SECTORS,  WHEN  THERE  ARE 

C 

C  ONLY  AAW  TYPE  THREATS. 

C 

C  SS .  VALUE  OF  THE  SURVIVAL  PROBABILITY  OF  THE  CARRIER  DUE 

C 

C  TO  ASW  THREATS  FROM  ALL  THE  ASW  SECTORS,  WHEN  THERE  ARE 

C 

C  ONLY  ASW  TYPE  THREATS. 

C 

C  T1 .  VALUE  OF  THE  SURVIVAL  PROBABILITY  OF  THE  CARRIER  DUE 

C 

C  TO  AAW  THREATS  FROM  ALL  THE  AAW  SECTORS,  WHEN  THERE  ARE 

C 

C  BOTH  AAW  AND  ASW  TYPE  THREATS. 

C 

C  T2 .  VALUE  OF  THE  SURVIVAL  PROBABILITY  OF  THE  CARRIER  DUE 

C 

C  TO  ASW  THREATS  FROM  ALL  THE  ASW  SECTORS,  WHEN  THERE  ARE 

C 
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BOTH  AAW  AND  ASW  TYPE  THREATS. 


SUBROUTINE  QDRTF(N,  X,  NF,  F,  UIP,  URP,  UFP) 

INTEGER  N ,NF,UIP (1) 

DOUBLE  PRECISION  X(N),  F,  URP(N,3) 

EXTERNAL  UFP 

DOUBLE  PRECISION  DEXP,  DSQRT 

DOUBLE  PRECISION  KA(I50) ,KS (150) ,  ASIGMA(150),  SSIGMA(150) 
REAL* 8  PI 

DOUBLE  PRECISION  AAWEXP(150),  AAWTOT(150),  AAWANS  U  "■) 
DOUBLE  PRECISION  ASWEXP(150),  ASWTOT(150),  ASWANS(150) 

REAL *8  RR,  SS,  Tl,  T2 

DOUBLE  PRECISION  AAWLL(150),  AAWUL(150) 

DOUBLE  PRECISION  ASWLL(150),  ASWUL(150) 

INTEGER  I,  IFAIL,  II 
INTEGER  C,  E,  G,  H,  M 
INTEGER  NSTOR(l) 

REAL *8  D01BAF 

DOUBLE  PRECISION  AAWFN1 ,  ASWFN1 

EXTERNAL  D01BAZ,  AAWFN1,  ASWFN1 

COMMON  /INDEX/  II 

COMMON  /V/  KA 

COMMON  /VV/  KS 

COMMON  /U/  ASIGMA 

COMMON  /W/  SSIGMA 

COMMON  /C/  AAWEXP 

COMMON  /D/  ASWEXP 

COMMON  /Y/  AAWLL,  AAWUL 

COMMON  /Z/  ASWLL ,  ASWUL 

COMMON  /ASECT/  G 

COMMON  /SSECT/  H 

COMMON  /SHIP/  M 

UIP(l)  =  NF 

NSTOR (1) =16 


DO  250  I  =  1,  N 

URP (1,2)  =  URP ( I , 1) *X ( I ) 
250  CONTINUE 
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*  *  *  DETERMINING  THE  THREAT  SECTORS  *  *  * 

G  =  0  :  NO  AAW  SECTORS 
H  =  0  :  NO  ASW  SECTORS 

X  X  x  x  x  * rxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

IF  THERE  ARE  NO  AAW  SECTORS  (I.E.  NO  AAW  THREATS  <  ,  THEN  G  =  C 
AND  WE  GO  TO  1111  TO  EVALUATE  FOR  ASW  THREATS. 

IF (G  .EQ.  0)  GO  TO  1111 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx*xx*x 

IF  THERE  ARE  NO  ASW  SECTORS  (I.E.  NO  ASW  THREATS),  H  =  0  AND 
WE  GO  TO  1112  TO  EVALUATE  FOR  AAW  THREATS. 

IF (H  .EQ.  0)  GO  TO  1112 

:XXXXXXXXXXXXXXX*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX-*:xXXXXXXX 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXxXXXXXXXXXXXXXXX 

IF  THERE  ARE  BOTH  AAW  AND  ASW  THREATS,  THEN  GO  TO  1113 

IF ( (G  .GT.  0)  .AND.  (H  .GT.  0))  GO  TO  1113 

’**X*XXXftXX***X**XXXXX*XXXXXXXXXXXftXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

GO  TO  1113 

'  »  CALL  INTEGRATION  ROUTINE  REPEATEDLY  OVER  G  AAW  SECTORS  *  *  * 
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1112  DO  65  II  =  1,G 
IFAIL=1 

AAWANS (II)  =  D01BAF (D01BAZ, AAWLL ( I I ) ,AAWUL(II) ,NSTOR(l) , 
1  AAWFN1 , IFAIL) 

65  CONTINUE 


C  *  *  CALCULATE  TOTAL  SURVIVAL  PROBABILITY  DUE  TO  ALL  A  AW  SECTORS 

C 

C 

RR  =  1.D+0 
DO  39  I  =  1,  G 

AAWTOT ( I )  =  AAWANS (I)  **  AAWEXP(I) 

RR  =  RR  *  AAWTOT (I) 

39  CONTINUE 
F  =  -RR 
GO  TO  999 


C  «  *  »  CALL  INTEGRATION  ROUTINE  REPEATEDLY  OVER  H  ASW  SECTORS  * 
C 

1111  DO  66  II  =  1 , H 
IFAIL  =  1 

ASWANS (II)  =  D01BAF (D01BAZ , ASWLL (II), ASWUL (II), NSTOR ( 1 ) , 

2  ASWFN1, IFAIL) 

66  CONTINUE 
C 
C 
C 

c 

c 

c  *  *  CALCULATE  TOTAL  SURVIVAL  PROBABILITY  DUE  TO  ALL  ASW  SECTORS  * 

C 

C 

SS  =  1.0D+0 
DO  38  I  =  1 , H 

ASWTOT ( I )  =  ASWANS  (I)  **  ASWEXP ( I ) 

SS  =  SS  *  ASWTOT ( I ) 
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38  CONTINUE 
F  =  -SS 
GO  TO  999 
C 
C 
C 
C 
C 

C  *  *  CALCULATE  SURVIVAL  PROBABILITY  DUE  TO  AAW  AND  ASW  SECTORS  *  * 

C 

C 

1113  DO  67  II  *  1 ,G 
IFAIL  *  1 

AAWANS (II)  =  D01BAF(D01BAZ,AAWLL (II) ,AAWUL(II) ,NSTOR(l) , 

1  AAWFN1, IFAIL) 

67  CONTINUE 


T1  «  1-D+O 
DO  34  I  *  If G 

AAWTOT(I)  *  AAWANS ( I )  **  AAWEXP(I) 
T1  *  T1  «  AAWTOT ( I ) 

34  CONTINUE 


DO  68  II  =  1 ,H 
IFAIL  »  1 

ASWANS (II)  *  D01BAF(D01BAZ,ASWLL(II) ,ASWUL(II) ,NSTOR(l) , 
2  ASWFN1, IFAIL) 

68  CONTINUE 


T2  *  1.D+0 
DO  35  I  *  If H 

ASWTOT(I)  »  ASWANS (I)  **  ASWEXP(I) 
T2  *  T2  *  ASWTOT ( I ) 

35  CONTINUE 


*  *  *  CALCULATE  OBJECTIVE  FUNCTION  *  *  * 

F  =  -(Tl  *  T2) 

GO  TO  999 
999  RETURN 
END 
C 
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AAWFN1 


C 
C 

c 
o 
c 
c 
c 

C  THE  FUNCTION  SUBPROGRAM  AAWFN1  CALCULATES  THE  SURVIVAL  PROBABILITY 
C 

C  OF  THE  CARRIER  IN  A  GIVEN  AAW  THREAT  SECTOR.  WE  FORMULATE  THE 
C 

C  INTEGRAND,  AAWFN1,  WHICH  IS  USED  BY  D01BAF. 

C 

C 


C... ****** ............... 

C 

C  PARAMETER  USAGE 

C 

c  a*******************************************  .*•«****•••••***•••*••* 

c  M 

C  AAWT . =  Jl(  1  -  PSHIP(I)  )  «  THE  PROBABILTIY  THAT  NO  SHIP 

C  rJ 

C  KILLS  THE  AAW  TARGET. 

C 

C  KLUGE ( I ) .  THE  OPTIMUM  INTERCEPT  DISTANCE  OF  SHIP  I  TO  A  TARGET. 

C 

C  PR (I) . THE  EFFECT  OF  THE  ITH  SHIP'S  DISTANCE  FROM  THE  CARRIER. 

C 

C  PSHIP(I) .  THE  PROBABILITY  THAT  SHIP  I  KILLS  A  TARGET. 

C 

C  R(I) . THE  DISTANCE  OF  THE  ITH  SHIP  FROM  THE  CARRIER,  I.E. 

C 

C  FROM  THE  ORIGIN. 

C 

C  THETA . THE  ANGLE  OF  APPROACH  OF  AN  AAW  TARGET.  IT  IS 

C 

C  UNIFORMLY  DISTRIBUTED  OVER  A  GIVEN  AAW  SECTOR. 

C 

C********** •••**•**•*••••**•• •••*•«•••**•**•• *•*•*•**»**.***••*•*••••*• 

c 


Appendix  C  The  Computer  Program  and  its  Documentation 


PAGE  192 


FUNCTION  AAWFN1 (THETA) 

DOUBLE  PRECISION  X(20) 

REAL *8  THETA 

REAL* 8  MINNUK 

REAL *8  SMALL 

REAL*8  AAWT 

INTEGER  C,  E,  1,  II,  M 

DOUBLE  PRECISION  KLUGE (150) ,  PSHIP(150) 

DOUBLE  PRECISION  AAVLL(150) ,  AAWUL(150) 

DOUBLE  PRECISION  KA(150),  ASIGMA(150) 

DOUBLE  PRECISION  R(150) ,  PR(150) 

DOUBLE  PRECISION  DSQRT ,  DEXP 

COMMON  /INDEX/  II 

COMMON  /V/  KA 

COMMON  /U/  ASIGMA 

COMMON  /SHIP/  M 

COMMON  /?/  AAWLL,  AAWUL 

COMMON  /X/  X 

COMMON  /MIN/  MINNUM 

COMMON  /SM/  SMALL 

*  *  *  COMPUTE  CARRIER  SURVIVAL  PROBABILITY  *  •  » 


AAWT  «  1.D+0 
DO  79  I  »  1,M 

C  ■  (2  •  I)  -  1 
E  ■  (2  •  I) 

R(I)  «  DSQRT ((X(C)*X(C))  +  (X(E)«X(E))) 

PR(I)  ■  i.D+0  -  DEXP(-KA(I)  *  R(I)) 

KLUGE ( I )  «  (X(C)*DSIN (THETA))  -  (X(E) *dCOS (THETA) ) 

IF  ( (-(KLUGE(I) •KLUGE(I)/ASIGMA(I) )  .LE.  MINNUM) 
l.OR.  (PR(I)  •  DEXP (-(KLUGE (I)  •  KLUGE(I)/ASIGMA(I) ) ) 

2.LE.  SMALL))  PSHIP(I)  -  0.D+0 

PSHIP(I)  »  PR(I)  •  DEXP (-(KLUGE (I)  *  KLUGE (I) ) /ASIGMA (I) ) 
AAWT  ■  AAWT  *  (1.D+0  -PSHIP(I)) 

79  CONTINUE 

*  •  •  THE  INTEGRAND  «  •  * 

AAWFN1*  ( (1.D+0)/ (AAWUL(II)-AAWLL(II) ) )  *  (1.D+0  -  AAWT) 
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J 

i 


ji 


< 


RETURN 

END 

C 

C 

c 

c 

c 

c**.«*. 

c 

C  ASUFN1 
C 

c*** ******* ********************* . . 

c 

c 

c 

C  THE  FUNCTION  SUBPROGRAM  AAVFN1  IS  USED  TO  CALCULATE  THE  SURVIVAL 
C 

C  PROBABILITY  OF  THE  CARRIER  IN  A  GIVEN  ASW  THREAT  SECTOR.  WE  FORMULATE 
C 

C  THE  INTEGRAND,  ASWFN1,  WHICH  IS  USED  BY  DOIBAF. 


C 

C 

C 

C 

C  PARAMETER  USAGE 

C 

C******************************»*«****************«*************»**»*»* 
C  M 

C  ASWT . »JI(  1  -  PSHIP(I)  )  »  THE  PROBABILITY  THAT  NO  SHIP 

C  '* 

C  KILLS  THE  ASW  TARGET. 

C 

C  KLUGE ( I ) . THE  OPTIMUM  INTERCEPT  DISTANCE  OF  SHIP  I  TO  A  TARGET. 

C 

C  PR(I) . THE  EFFECT  OF  THE  ITH  SHIP'S  DISTANCE  TO  THE  CARRIER. 

C 

C  PSHIP(I) .  THE  PROBABILITY  THAT  SHIP  I  KILLS  AN  ASW  TARGET. 

C 

C  R(I) .  THE  DISTANCE  OF  THE  ITH  SHIP  FROM  THE  CARRIER,  I.E. 

C 

C  FROM  THE  ORIGIN. 

C 

C  THETA . THE  ANGLE  OF  APPROACH  OF  AN  ASW  TARGET.  IT  IS  UNIFORMLY 

C 

C  DISTRIBUTED  OVER  A  GIVEN  ASW  SECTOR. 
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FUNCTION  ASWFN1 (THETA) 

DOUBLE  PRECISION  X(20) 

REAL *8  THETA 

REAL*8  HINNUX 

REAL *8  SMALL 

REAL* 8  ASWT 

INTEGER  C,  E,  I,  II,  M 

DOUBLE  PRECISION  KLUGE(150),  PSHIP(150) 

DOUBLE  PRECISION  ASWLL(l50),  ASWUL(150) 

DOUBLE  PRECISION  KS(150),  SSIGMA(150),  R(150) ,  PR(150) 

DOUBLE  PRECISION  DSQRT ,  DEXP 

COMMON  /INDEX/  II 

COMMON  /VV/  KS 

COMMON  /V/  SSIGMA 

COMMON  /SHIP/  M 

COMMON  /Z/  ASWLL,  ASWUL 

COMMON  /X/  X 

COMMON  /MIN/  MINNUM 

COMMON  /SM/  SMALL 


*  *  »  COMPUTE  CARRIER  SURVIVAL  PROBABILITY  *  *  * 

ASWT  ■  l.D+O 
DO  81  I  »  1>M 

C  »  (2  »  I)  -  1 
E  *  (2  *  I) 

R(I)  «  DSQRT  ((X(C)*X(C))  +  (X(E)*X(E))) 

PR(I)  »  l.D+O  -  DEXP(-KSU)  •  R(I) ) 

KLUGE ( I )  «  (X(C)«DSIN (THETA))  -  (X(E) *DCOS (THETA) ) 

IF  ( (-(KLUGE(I) *KLUGE(I)/SSIGMA (I) )  .LE.  MINNUM) 
l.OR.  (PR(I)  *  DEXP (-(KLUGE (I)  •  KLUGE (I) /SSIGMA (I) ) ) 

2.LE.  SMALL))  PSHIP(I)  »  0.D+0 

PSHIP(I)  «  PR(I)  •  DEXP (-(KLUGE (I)  «  KLUGE (I) )/SSIGMA(I) ) 
ASWT  »  ASWT  *  (l.D+O  -PSHIP(I) ) 

81  CONTINUE 


•  «  *  THE  INTEGRAND  *  *  • 

ASWFN1*  ( (l.D+O)/ (ASWUL (II) -ASWLL (II) ) )  *  (l.D+O  -  ASWT) 


Appendix  C  Th«  Computer  Program  and  its  Documentation 


PAGE  195 


RETURN 

END 


•V 
»•  * 


•**  SUBROUTINE  READ3  READS  THE  VALUES  FROM  THE  DATA  FILE  *** 


READ3 


SUBROUTINE  READ3 (G,H,K,NN ,D,X,KA,KS,ASIGKA,SSIGHA,AAWLL, AAWUL, 
1AAWEXP f ASWLL , ASVUL , ASWEXP , HINNUH , SMALL ) 

INTEGER  G,H,M,NN 
REAL'S  MINNUM , SMALL 

DOUBLE  PRECISION  D(20) ,  X(20),  ASIGMA(150) ,  SSIGMA(150) 

DOUBLE  PRECISION  AAWLL(150),  AAWUL (150) ,  ASWLL ( 150) ,  ASWUL(150) 
DOUBLE  PRECISION  AAWEXP(150) ,  ASWEXP (150) 

DOUBLE  PRECISION  KA(150) ,  KS(150) 


*  •  «  READ  the  IMPLEMENTATION  DEPENDENT  CONSTANTS  * 


READ(82,l)  MINNUM,  SMALL 
1  FORMAT (/D13. 6, 2X,E7.1) 


«  *  *  READ  THE  TOTAL  NUMBER  OF  COORDINATES  »  *  » 


READ (82, 2)  NN 
2  FORMAT (///I2) 
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«  *  *  READ  THE  SCALE  ARRAY 


*  *  * 


DO  10  I  *  1,NN 
READ (82,3)  D ( I ) 

3  FORMAT (D9. 3) 

10  CONTINUE 
C 

C 

C  *  *  »  READ  THE  INITIAL  CONDITIONS  «  *  « 

C 

C 

DO  11  I  »  1,NN 
READ (82 ,4)  X(I) 

4  FORMAT (D9. 3) 

11  CONTINUE 
C 

C 

C  *  •  «  READ  THE  TOTAL  NUMBER  OF  SHIPS  *  *  * 

C 

C 

READ (82 .5)  M 

5  FORMAT (///I2) 

C 

C 

c  »  «  •  READ  THE  AAV  DISTANCE  SCALE  FACTORS  *  «  » 

C 

C 

DO  12  I  «  1*  M 

READ (62, 6)  KA(I) 

6  FORMAT (D9. 3) 

12  CONTINUE 
C 

C 

C  •  •  «  read  THE  ASW  DISTANCE  SCALE  FACTORS  *  *  * 

C 

C 

DO  121  I  -  1,M 

READ (82, 61)  KS(I) 

61  FORMAT (D9. 3) 

121  CONTINUE 
C 
C 

C  *  *  *  read  THE  AAV  RANGE  FACTORS  *  *  * 

C 

C 
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DO  13  I  *  1,M 

READ (82, 7)  ASIGKA(I) 
7  FORM AT (D9. 3) 

13  CONTINUE 


*  *  *  READ  THE  ASW  RANGE  FACTORS  «  *  * 


DO  14  I  *  1»M 

READ (82 ,8)  SSIGMA(I) 
8  FORMAT (D9. 3) 

14  CONTINUE 


*  «  «  READ  THE  TOTAL  NUMBER  OF  AAW  THREAT  SECTORS  *  *  * 


READ (82, 9)  G 
9  FORMAT (///I2) 


IF(G  .EQ.  0)  GO  TO  555 


•  *  •  READ  THE  LOWER  LIMITS  OF  ALL  THE  AAW  SECTORS  *  *  * 


DO  15  I  »  1,6 

READ (82, 100)  AAWLL(I) 
100  FORMAT(D9.3) 

15  CONTINUE 


»  •  *  READ  THE  UPPER  LIMITS  OF  ALL  THE  AAW  SECTORS  *  *  * 


DO  16  I  >  1 ,  G 

READ (82, 101)  AAWUL(I) 
101  FORMAT (09. 3) 

16  CONTINUE 


•  •  •  READ  THE  NUMBER  OF  TARGETS  IN  EACH  AAW  SECTOR  *  *  * 
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C 

c 

DO  17  I  ■  1 ,G 

READ(82,102)  AAWEXP(l) 

102  FORMAT (D9. 3) 

17  CONTINUE 

C 

C 

c  *  *  *  read  the  total  number  of  asw  threat  sectors  *** 

c 

c 

555  READ (82, 103)  H 

103  FORMAT (///I2) 

C 

C 


IF(H  .EQ.  0)  GO  TO  999 
C 
C 

C  *  *  *  READ  THE  LOWER  LIMITS  OF  ALL  THE  ASW  SECTORS  «  *  * 

C 

C 

DO  18  I  ■  1,H 

READ (82, 104)  ASWLL(I) 

104  FORMAT (D9. 3) 

18  CONTINUE 
C 

C 

C  *  *  «  READ  THE  UPPER  LIMITS  OF  ALL  THE  ASW  SECTORS  «  *  * 

C 

C 

DO  19  I  «  1,H 

READ (82, 105)  ASWUL(I) 

105  FORMAT (D9. 3) 

19  CONTINUE 
C 

C 

C  •  •  *  READ  THE  NUMBER  OF  TARGETS  IN  EACH  ASW  SECTOR  *  *  * 

C 

C 

DO  20  I  ■  1,H 

READ (82, 106)  ASWEXP ( I ) 

106  FORMAT (D9. 3) 

20  CONTINUE 
C 
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D.l  SAMPLE  OUTPUT 


Sample  output  of  a  specific  experiment  follows.  Included  is  the  associated 
data  file,  KSH  DT. 
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MINNUM  SMALL 

-0.174673D+03  1.0E-75 

C 

*  *  *  #  OF  COORDINATES  :  NN  «  *  * 


C 

4 


0.100D+01 

=  D(l) 

0.100D+01 

=  D  (2) 

0.100D+01 

=  D(3) 

0 • 100D+01 

=  D(4) 

0.170D+02 

=  X(l) 

0.170D+01 

=  X  (2) 

0.100D+02 

=  X(3) 

0.380D+02 

=  X  (4) 

C 

«  *  *  TOTAL  # 

OF  SHIPS  :  M 

C 

2 

0 . 200D+00 

-  KA(1) 

0.200D+00 

=  KA(2) 

0.100D+00 

=  KS  (1) 

0.100D+00 

=  KS  (2) 

0.100D+04 

=  ASIGMA(l) 
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0.850D+03  =  ASIGMA (2) 

0.500D+03  =  SSIGMA(l) 

0.800D+04  *  SSIGMA (2) 

C 

*  *  *  #  OF  AAW  SECTORS  :  G  *  *  « 

C 

1  #  OF  AAW  SECTORS 

O.OOOD+OO  «  AAWLL 

0.600D+02  *  AAWUL 

0.400D+01  *  AAW EXP 

C 

*  *  *  #  OF  ASW  SECTORS  :  H  «  »  « 

C 

1  #  OF  ASW  SECTORS 

0 . 300D+02  *  ASWLL 

0.900D+02  *  ASWUL 

0.100D+01  a  ASWEXP 
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